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SCATTERING FACTORS FOR SOME OF THE HEAVIER ATOMS! 


BEATRICE H. WorsLEY 


ABSTRACT 


A program for calculating X-ray atomic scattering factors from the radial 
wave functions has been written for the IBM 650 installation at the University 
of Toronto. It has been applied to the results of self-consistent field calculations 
previously performed at this University on the FERUT computer. Results are 
given for Ne, V**, Kr, Agt, and Pb*+**. The results for Ne and V++ are com- 
pared with those calculated by Freeman using Allen’s wave functions for Ne 
and Hartree’s approximate wave functions for Vt*. 


INTRODUCTION 


As new Hartree-Fock radial wave functions become available, it becomes 
possible to derive functions which can more readily be compared with ex- 
perimental results. Such a function is the X-ray scattering factor for coherent 
radiation, which is of particular interest to crystallographers. It has been 
described by Hoerni and Ibers (1954). This function has been calculated for 
several specific cases by hand and also with the aid of automatic digital 
computers. Just recently it has been calculated for 45 different atoms and 
ions by Freeman (1959) using the Whirlwind computer at M.I.T. 

A program has been written at the University of Toronto for the IBM 650 
installation to calculate the atomic scattering factor for a spherical charge 
distribution, given the appropriate set of radial wave functions. This program 
has been applied to the self-consistent fields with exchange for Ne (Worsley 
1958a) and for V+*, Kr, and Agt (Worsley 1958). It has also been applied 
to the less accurate self-consistent field without exchange for Pb++* (Hart and 
Worsley 1959). 

THE PROGRAM 


Assuming that the electron density is spherically symmetric and averaging 
over all orientations, Hoerni and Ibers (1954) give, for the definition of the 
atomic scattering factor: 


f(s) = f° > {P*(r) sin sr/sr}dr, 


0 all 
electrons 


‘Manuscript received April 30, 1959. j : 
Contribution from the Computation Centre, McLennan Laboratory, University of 
Toronto, Toronto, Ontario. 
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where the P(r)’s are the radial wave functions and s = 4x\~' sin @ is the 
magnitude of the radiation vector in reciprocal space, s being treated as a 
parameter. 

To be readily applicable to the wave functions calculated by FERUT, the 
present program was designed to evaluate this function in the form 


1005p 


f(s) = a {P?(p) sin sr/s}dp, 


p=0 a! 
electrons 


where the P’s are supplied at equal intervals of the logarithmic variable 
p = Incr. For each application it should be verified that p = 0 is an adequate 
approximation to p = —~, and p = 100ép to p = + in the sense that 
contributions to the integral outside the range 0 < p < 1006p are negligible. 
Provision is made for introducing ¢ and 6p as well as s into the program as 
parameters. Furthermore, the program will handle a configuration with as 
many as 16 different wave functions. Newton’s three-eighth’s rule is used to 
evaluate the integral. It is modified in such a way that, in effect, the half- 
interval is used, the values of the integrand at the mid-points being evaluated 
by linear interpolation. 


RESULTS 


Results of the current application of this program are given in Table I. 
The errors are indicated, on the assumption that the input wave functions 


TABLE I 


Values of the atomic scattering factors 





Agt, Pbt++,* 





Ne, Vit, Kr, 

sin 6/x A> f+0.003 f+0.01 f+0.01 f+0.0 f+0.1 
0.00 10.000 21.00 36.00 46.00 79.0 
0.05 9.832 20.70 35.30 45.38 78.3 
0.10 9.351 19.86 33.45 43.68 76.0 
0.15 8.643 18.62 31.00 41.20 72.6 
0.20 7.804 17.15 28.47 38.31 68.8 
0.25 6.926 15.60 26.15 35.33 64.8 
0.30 6.077 14.10 24.12 32.48 60.9 
0.35 5.300 12.72 22.37 29.89 57.2 
0.40 4.614 11.51 20.83 27.63 53.8 
0.50 3.531 9.61 18.15 24.02 47.9 
0.60 2.786 8.33 15.81 21.42 43.1 
0.70 2.292 7.47 13.73 19.48 39.0 
0.80 1.968 6.86 11.95 17.89 35.5 
0.90 1.754 6.38 10.46 16.48 32.3 
1.00 1.607 5.94 9.26 15.14 29.4 
1.10 1.500 5.52 8.32 13.86 26.7 
1.20 1.416 5.10 7.59 12.63 24.4 
1.30 1.344 4.69 7.03 11.48 22.3 
1.50 1.216 3.91 6.25 9.50 19.1 
1.70 1.099 3.24 5.69 7.98 16.8 
1.90 0.981 2.69 5.23 6.90 15.2 





*Wave functions without exchange. 
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are correct. It should be noted that the sine term modulates somewhat rapidly 
for the higher values of s. Nevertheless, it seemed worth while to include 
these in the calculations. 

In all cases of functions derived from atomic wave functions, the accuracy 
of the latter should be taken into account. This accuracy depends upon the 
closeness of the formulation to the physics of the system and upon the numerical 
errors introduced in the course of the calculation. In the case of the results 
in Table I, the errors of calculation of the wave functions are overridden by 
the errors in the calculation of the scattering factors, so that the results 
should reflect the closeness of the Hartree-Fock formulation to the physical 
picture. 

The values for Ne and V+* have been compared with those recently pub- 
lished by Freeman (1959). In the case of Ne, Freeman used L. C. Allen’s un- 
published wave functions, and the maximum discrepancy in the scattering 
factors is 0.4%. For V++, Freeman used the approximate wave functions 
published by Hartree (1956), and the discrepancies are correspondingly 
greater. These amount to as much as 7% at sin 6/A = 1.00, or as much as 
0.80 A-! units at sin 0/A = 0.40. The present values reflect the higher accuracy 
of the FERUT wave functions. Finally, it has been reported by Ibers (1959) 
that the results for Pb+++ are in good agreement with values based on cal- 
culations for the Thomas—Fermi-—Dirac model. 
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THE ONSET TIMES OF FORBUSH-TYPE COSMIC RAY 
INTENSITY DECREASES! 


A. G. FENTOoN,? K. G. McCrackKeEn,? D. C. Rose, AND B. G. WILSON? 


ABSTRACT 


The onset times of a number of Forbush-type decreases observed at four widely 
spaced stations are compared, and it is shown that appreciable differences occur. 
The stations selected were Hobart, Mawson, Ottawa, and Sulphur Mountain. 
It was found that a consistent pattern is obtained for the events studied when 
the onset times are plotted as a function of the direction of maximum sensitivity 
of the recorders relative to the earth-sun line. This is interpreted as being due 
to a directional anisotropy that exists in the mechanism producing the decreases, 
at least in the early stages. The depression occurs first for particles arriving from 
directions between 30° and 120° west of the earth-sun line. The relation between 
these observations and geomagnetic disturbances and the quiet-day daily varia- 
tion is discussed. 


INTRODUCTION 

It was mentioned briefly in a recent paper (Fenton, McCracken, Rose, and 
Wilson 1959) that an examination of the hourly intensities at widely separated 
stations has revealed obvious differences in response to Forbush-type de- 
creases, and that these differences appear to be at least partly due to primary 
anisotropies which are detected at different times and directions in space as 
the earth rotates. In this report we shall consider in some detail the time 
dependency of the commencements of a number of events of this type ob- 
served at the stations operated by the Ottawa and Hobart groups during 
recent years of high solar activity. During this period data were obtained 
from high counting rate recorders developed for the International Geophysical 
Year 1957-58, enabling the analysis of individual events at reasonable counting 
statistics, thus avoiding recourse to the procedure of averaging over a number 
of events. 

Data from four widely spaced stations, namely, Hobart, Mawson, Ottawa, 
and Sulphur Mountain, will be presented and examined for evidence of a 
consistent pattern in the intensity-time changes at the commencement of 
Forbush-type decreases. If it exists, such a pattern would be expected to 
provide important information concerning the space-time behavior of the 
agency responsible for these transient cosmic ray intensity decreases. At the 
present state of knowledge of the mechanisms involved we are of the opinion 
that much can be gained by a phenomenological approach to these problems, 
using the currently available data of adequate statistical accuracy in hourly 
values, in order to establish firmly the essential features upon which a sound 
physical explanation may be based. 

'1Manuscript received June 2, 1959. 

Contribution from the Department of Physics, University of Tasmania, Hobart, and the 
Division of Pure Physics, National Research Council, Ottawa, Canada. 

Issued as N.R.C. No. 5274. 


2Department of Physics, University of Tasmania, Hobart, Tasmania. 
3National Research Council Sulphur Mountain Laboratory, Banff, Alberta. 
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THE DEPENDENCE OF ONSET TIMES UPON STATION CO-ORDINATES 

Examination of the hourly intensities from widely separated stations, 
plotted as a function of universal time, reveals that there are obvious differences 
in the observed onset time of many Forbush-type decreases. Also, these 
differences seem to depend more on the station longitudes than on their 
latitudes. Sample events to illustrate this are shown in Fig. 1, in which are 
plotted the intensities observed at a number of stations. The co-ordinates of 
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Fic. 1. The pressure-corrected data obtained during the initial phases of the Forbush-type 
decreases occurring on 21 October, 1957, and 26 November, 1957. 
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stations* operated by the Ottawa and Hobart cosmic ray groups have been 
set out previously (Fenton, McCracken, Rose, and Wilson 1959). It may be 
noted that there is a close similarity between intensity-time plots for data 
from the stations lying within a narrow longitude zone, namely, Hobart, 
Macquarie Island, and Lae. On the other hand, the Mawson, Ottawa, and 
Sulphur Mountain results have quite a different appearance. 

Associated with some events there is a clearly distinguishable short-lived 
decrease which occurs before the main decrease. Cases. of this phenomenon 
are included in Fig. 1; for example, it occurred at Hobart during the October 21 
event and at Mawson during the November 26 event. One of these short- 
lived transient decreases has been described recently (McCracken and Parsons 
1958) as a “‘preliminary decrease’’ associated. with the October 21 event. As 
will be pointed out below, this type of intensity change has been observed on 
a number of occasions at one or more of our stations, so it is by no means a 
rare phenomenon. In fact, it seems likely that the mechanism responsible for 
this type of variation is present, to a greater or lesser extent, during the early 
stages of all Forbush-type events. 

Following the technique used previously (Fenton, McCracken, Rose, and 
Wilson 1959) to study the superposition of transient decreases, the event of 
November 26, 1957, has been examined as it appears in the average of the 
intensities recorded at Hobart, Mawson, Ottawa, and Sulphur Mountain. 
This procedure is expected to give an approximation to the decrease that 
would be observed by an omnidirectional (47) detector placed at the position 
of the earth, and thus able to detect temporal changes in the depressive 
mechanism, independently of scanning effects. In the averaging process each 
station was allotted equal weight, irrespective of counting rate. No allowance 
was made for the effect of the different station altitudes upon the observed 
event amplitude. Although recent work has shown that correction for altitude 
is desirable for some studies (McCracken and Johns 1959), neglect of this is 
not serious in the present analysis. The average response at the four stations 
during the November 26 event is plotted in Fig. 2, which shows that the 
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Fic. 2. Average nucleon intensities observed at four stations during November 26 event, 
1957, expressed as percentages of the mean intensity observed on November 21. 


*Co-ordinates of stations reported here: 
Latitude Longitude 


Ottawa 45.4N. 75.6W. 
Sulphur Mountain 51.1N. 115.6 W. 
Hobart 42.8S. 147.5E. 
Mawson 67.6S. 62.9 E. 
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event may be regarded as a rather slow change upon which is superimposed a 
more rapid short-lived event. This type of analysis has been made for other 
events, and generally leads to a smoothing of the intensity-time plot. This 
suggests that the variation in the time of commencement, the rate of change 
of intensity, and the general appearance of the intensity—time plot for Forbush- 
type decreases are due to the directional sensitivity of the recorders and 
asymmetries in the agency responsible for the decreases, made evident by the 
rotation of the earth. In other words, the integral effect at the earth due to 
the cosmic ray intensity-reducing agency varies more smoothly with time 
than is apparent from the records at a single station. 

Intensity-time plots similar to those of Fig. 1 have been prepared for 
events observed between August 1956 and December 1958. Information from 
events of amplitude about 4% or greater is summarized in Table I. During 


TABLE I 


The times of commencement of Forbush-type decreases (or of the preliminary phases of 
these events when identified) are given for four stations, together with the magnetic storm 
sudden commencements to which the decreases may be related 


Hour of commencement, U.T. 


Sulphur = S.S.C., Preliminary decrease 


Date of event Hobart Ottawa Mawson Mt. U.T. observed at: 
20 Aug. 1956 16 09 (T) 01 (21st) Hobart, Mawson, 
(T) Ottawa 
9 Nov. 1956 02(10th) 23 22:(T) 2000 
25 Dec. 1956 11 03 12(T) 0800 
21 Jan. 1957 21 21 19 (T) 1230 
10 Mar. 1957 09 05 00 (T) 0030 
17 Apr. 1957 08 00 ll 1130 
22 June 1957 08 04 04 07 
3 Aug. 1957 = 13 24 16 15 1600 ‘Four stations 
29 Aug. 1957 21 23 20 21 1930 
29 Sept. 1957 04 06 (28th) 02 04 0030 Mawson, Ottawa 
21 Oct. 1957 19 09 23 02 (22nd) 2330 Hobart, Ottawa 
26 Nov. 1957 07 05 04 06 Mawson 
11 Feb. 1958 03 02 03 02 0130 
25 Mar. 1958 12 16 14 19 1630 
9 May 1958 10 08 05 09 
29 May 1958 12(28th) 04 07 04 Hobart 
28 June 1958 14 09 15 01 
8 July 1958 12(7th) 10(7th) 21(7th) 10 0800 Hobart, Ottawa, 
Mawson 
17 Aug. 1958 06 05 06 01 0730 Hobart, Ottawa, 
Sulphur Mt. 
24 Aug. 1958 19(23rd) 04 01 03 0300 Hobart 
15 Sept. 1958 16 09 23 03 (16th) Ottawa 
22 Oct. 1958 24 23 03 (23rd) 23 0330 


Note: Times corresponding to neutron monitor data are given except where indicated by the letter T for meson 
telescope. Blank spaces indicate equipment not operating, S.S.C. not observed or preliminary decrease not identified. 
Hour 24 refers to the intensity observed between 2300 and 2400 hours U.T. Hour 00 indicates hour 24 of the 


previous day. 


preparation of this table the onset time of a Forbush-type decrease has been 
defined as the first hour of the decreasing phase of the event during which 
the intensity was more than two standard deviations below the pre-event 
level. The onset time listed is that of the preliminary decrease where such a 
decrease is clearly connected with the event. It is sometimes difficult to decide 
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upon a pre-event level if there are concurrent intensity changes other than 
statistical fluctuations. In such cases a smooth curve was drawn through the 
points prior to the first obviously decreasing phase of the event. On some 
occasions there is a pronounced maximum intensity just before the decrease 
begins. It is not known whether this is a significant feature or whether it 
arises simply because of recovery from earlier transient decreases. Through- 
out the present investigation nucleon intensity data corrected for pressure 
have been used in preference to meson data, although the latter have been 
employed to provide supporting evidence when determining onset times, or 
to permit the use of events that occurred before neutron monitors were in 
operation at all four stations. There is an indication that the onset times at 
a given station are not always the same for the nucleon and meson components, 
but the time difference is small, if genuine, and is not likely to influence the 
conclusions of this paper. 

Onset-time differences amounting to many hours are apparent from Table | 
and Fig. 1. Errors in the determination of onset time have not been assigned, 
but would appear to vary from 0 to 2 or 3 hours depending on the abruptness 
of the event and the counting rate. 

THE DIRECTIONAL SENSITIVITY OF THE RECORDERS 

In what follows, we take the view that the asymmetries resulting in the 
rapid changes in intensity, as the earth rotates, are situated outside the region 
occupied by the earth’s magnetic field. The justification for this assumption 
is the observation of a number of preliminary depressions and Forbush de- 
creases beginning prior to the commencement of the accompanying magnetic 
storm. This suggests that the corpuscular matter responsible for the magnetic 
storm, while outside the region occupied by the geomagnetic field, is able 
to affect the cosmic ray intensity at the earth. 

The particles detected by a neutron monitor are due to primaries arriving 
from within a cone of acceptance of rather small solid angle compared with 
4x, so that as the earth rotates the recorder scans a region of space determined 
by the location of the recording station and geophysical factors. Thus a daily 
survey of the intensity arriving from different directions relative to the earth— 
sun line is obtained by selecting data from stations at different latitudes. 
If sufficient stations at different longitudes as well as latitudes are selected, 
a practically continuous survey of all regions is possible. 

When determining the cone of acceptance of a recorder it is necessary to 
take into account the shape of the primary spectrum and the deflection of 
charged particles in the earth’s magnetic field. In addition there are effects 
due to the properties of the nucleon cascades in the atmosphere (the specific 
yield function) and geometrical factors. 

Taking all these factors into account, the dependence of counting rate upon 
direction has been determined for a number of stations, the dependence of 
counting rate upon geographic longitude for the stations considered in this 
paper being given in Fig. 3. In the calculation of these curves, the published 
data on the deflection of cosmic rays in a dipole field (Brunberg 1953; Brun- 
berg and Dattner 1953; Jory 1956), and the coupling constants appropriate 
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to the neutron component (Dorman 1957), have been used. A description of 
the derivation of these curves will be given in a later paper. 

It will be seen from Fig. 3 that the direction of maximum response of a 
neutron monitor is not at a fixed angular distance from the meridian plane 
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THE ORDINATE SCALES ARE ARBITRARY 





60° 40° 20° 0° 20° 40° 60° 80° 100° 
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LONGITUDE RELATIVE TO OBSERVATORY 


Fic. 3. The contributions made to the total counting rates of four different neutron moni- 
tors as functions of asymptotic longitude. Longitude is measured relative to the meridians 
through the monitors. The ordinate scales are arbitrary. 


through the station. The longitude of the directions of maximum response 
relative to the station meridian is set out in Table II for the four stations of 


TABLE II 


Geographic longitudes of observatories and directions 
of maximum response measured east from the 
Greenwich meridian 


Direction of 
ae | maximum 
Station | Observatory response 


~_ Longitudinal | 





Mawson 63° E. 68° E. 
Hobart 147 207 
Sulphur Mt. 245 270 
Ottawa 284 319 
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interest. It is clear that a region in space from which the intensity is depressed 
will be scanned at different local times by different stations. For this reason 
we prefer to express the results presented in this paper in universal time 
rather than local time, and to relate the observations to directions in space 
as described below. 
THE TIME AND DIRECTIONAL DEPENDENCE OF THE FORBUSH-TYPE 
DECREASE MECHANISM 


Three well-defined Forbush events are represented in Fig. 4 as follows: 


TO SUN 
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Fic.4, (A) Showing the directions of greatest response of four neutron monitors at 0000 U.T. 
The positions of the recorders themselves are indicated by the solid circles. The convention 
whereby directions are designated as being west or east of the earth-sun line is illustrated. 

(B, C, D) Surveys of the directional and temporal dependence of intensity during three 
Forbush-type decreases. The letters indicate from which recorder the adjacent data were 
obtained. Thus M—Mawson, H—Hobart, S—Sulphur Mountain, O—Ottawa, C—Climax. 
The times at which magnetic storm sudden commencements occurred are indicated, 
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In these diagrams the ordinates represent the progress of time and the 
abscissae denote asymptotic directions east or west of the earth-sun line (see 
Fig. 4A and Table II for a definition of these directions). As the intensities 
are sampled hour by hour while the earth rotates a small square representing 
the direction of maximum sensitivity of a given recorder moves across the 
diagram from the lower left to upper right. An open square indicates no signi- 
ficant depression below the pre-event level, while a blackened square identifies 
times and directions for which the intensity was depressed below this level 
by more than two standard deviations. The contour lines serve to distinguish 
depressed from non-depressed zones of the diagrams. 

It is clear that on the three occasions represented in Fig. 4 the onset be- 
havior of the events had similar features. The most striking feature is that 
the radiation arriving from directions out to about 120° west of the earth-sun 
line is depressed first, while that from 0 to 120° east is depressed last. 


IDEALIZED TIME-DIRECTION DIAGRAM FOR FORBUSH-TYPE DECREASES 


Two of the events depicted in Fig. 4 have been examined in greater detail. 
These occurred on October 21 and November 26, 1957, and were well-defined 
cases of variable onset behavior from station to station. Instead of the 
“‘depressed’’ or ‘‘not depressed’’ representation we shall now indicate the 
extent of depression below the pre-event level. This is illustrated in Fig. 5 in 
which different symbols are used for each 2% interval between 2% and 12% 
depressed intensity. To allow for the altitude dependence (McCracken and 
Johns 1959) the Sulphur Mountain data have been multiplied by the factor 
0.8 to make them comparable with those from the other stations. Figure 5 
points to a systematic behavior for these two events, the values of greatest 
depression avoiding the directions between 0 and 90° east. The intensity due 
to particles arriving from east was the last to be depressed and it suffered the 
least depression during the events. 

On the basis of these two well-defined events the idealized time-direction 
diagram shown in Fig. 5C was prepared. From this diagram we may predict 
the following features to be expected on other occasions: 

(a) When they occur, short-lived preliminary depressions will be observed 
at times characteristic of the station. According to the diagram these times 
are centered about: Hobart 1800-2000 U.T., Mawson 0400-0600 U.T., 
Ottawa 1100-1300 U.T., Sulphur Mountain 1400-1600 U.T. The duration 
of the preliminary depression will be determined to an extent by the shape of 
the directional sensitivity curve appropriate to the station. 

(b) Intensity maxima occurring during the progress of the decrease will 
be observed at characteristic times. These are: Hobart 2300-0400 U.T., Maw- 
son 0800-1300 U.T., Ottawa 1500-2000 U.T., Sulphur Mountain 1900-2400 
U.T. 

In Table III these predicted features are compared with the observations 
for five well-defined but less striking events than those used in Fig. 4. The 
good agreement leads us to conclude that the idealized time—direction diagram 
adequately describes the structure of at least eight of the Forbush-type 
decreases examined. 
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Fic. 5. (A, B) Detailed surveys of two Forbush-type decreases. The symbols are: cross = 2 
to 4% intensity depression; open circle = 4 to 6% depression; small solid circle = 6 to 8% 
depression; large solid circle = 8 to 10% depression; solid square = 10 to 12% depression. 

(C) Survey of the idealized Forbush-type decrease, showing the manner in which the depres- 
sion of intensity depends upon direction and time. The directions of greatest depression are 
indicated by dense dotting. 

(D) Survey of two great Forbush-type decreases. Hy corresponds to the Hobart meson 
data; Hy, the Hobart neutron data. Times of sudden storm commencement are indicated. 


Two great transient decreases occurred during 1957, namely, those beginning 
on January 21 and August 29. The nucleon intensity decrease exceeded 10% on 
these occasions, but the observed onset times at our stations were not very 
different, and there were no preliminary decreases. Time-direction plots for 
these two events are shown in Fig. 5D. As the Sulphur Mountain and Mawson 
neutron monitors were not in operation until after the January event other 
data have been used. Although there were evidently no strongly developed 
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TABLE III 


The times of occurrence of prominent features of Forbush-type decreases, compared with 
those predicted from the idealized model 


~~ Feature 


























Short-lived (preliminary) decrease Intensity maximum observed during 
Event prior to main decrease fluctuations after onset 
Date of Sulphur Sulphur 
commencement Hobart Mawson Ottawa Mt. MHobart Mawson Ottawa Mt. 
20 Aug. 1956 17-20 01-06 09-11 N.D. 00-05 07-16 16-20 N.O. 
(Meson) (Meson) (Meson) 
10 Mar. 1957 ND: N.D. N.D. N.O. 20-01 05-12 14-17 N.O. 
3 Aug. 1957 18-22 03-10 08-10 16-20 20-04 11-14 N.D. 05-07 
29 Sept. 1957 N.D. 04-06 09-11 N.D. 23-04 08-14 09-15 22-24 
8 July 1958 15-18 01-05 10:14 ND. ND. ND. ND. ND. 
17 Aug. 1958 16-20 N.D. 09-13 13-17 23-02 N.D. 20-22 20-33 
Predicted central 
15-20 19-24 


values 18-20 04-06 11-13 14-16 23-04 08-13 








Note: Tabulated values are in universal time. Meson telescope data are used where nucleon data are not avail- 
able. N.D. indicates that the feature was not well defined in the data, and N.O. that the station was not operating 
at this time. 


anisotropies present at the time, these events appear to be consistent with the 
idealized diagram, since the intensity was depressed in the westerly before the 
easterly directions. 


DISTRIBUTION OF ONSET TIMES OF FORBUSH-TYPE DECREASES 

The results given in the previous section indicate that at least on a number 
of occasions the cosmic rays approaching the earth during the early stages 
of a Forbush-type event are reduced in intensity in an asymmetrical fashion. 
The suggestion from these results is that the intensity-depressing agency is 
located at first in a direction between 0-120° west of the sun-earth line, and 
that it later extends to influence the cosmic ray flux arriving from all direc- 
tions. It is obviously desirable to determine whether this form of behavior is 
common to all events. Evidence strongly in support of this is presented in 
Fig. 6. In this figure the earliest onset times listed in Table I for each event are 





000 900 90 160 
WEST EAST 


Fic. 6. The earliest time of commencement for the events listed in Table I is indicated by 
a black square in the time-direction diagram. When the commencement times at two stations 
were the same both have been shown. The number 2 alongside a square indicates that two 
events had the same commencement time. 
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plotted on a time-direction diagram. Only those events for which data were 
available from the four stations have been used. It is remarkable that although 
(as would be expected) the onset times are distributed throughout the 24 
hours, the directions in which the events were first observed to begin are by 
no means evenly distributed. Of the 22 onset times included in the figure none 
corresponds to directions from 45° west to 60° east of the sun-earth line. 

We conclude from the above analysis of 22 events that at least during the 
present cycle of high solar activity the Forbush-type decrease mechanism has 
associated with it a consistent asymmetrical pattern about the earth-sun 
line. One may summarize as follows: 

(a) The onset times are not simultaneous at all points on the earth. 

(b) The intensity between about 30° and 120° west of the earth-sun line 
is the first to be depressed. 

(c) The intensity is depressed last from directions between about 0° and 90° 


east. 
(d) The maximum depression occurs about 90° west. 


CORRELATION WITH GEOMAGNETIC DISTURBANCE 


Considerable attention has been devoted in the past to the connection be- 
tween geomagnetic activity and cosmic ray intensity variations. Usually the 
counting statistics have been inadequate to justify the study of individual 
events; consequently superposed epoch methods have been used. Conflicting 
results have been reported concerning the correlation between the two pheno- 
mena, but most investigators agree that the relationship is not that of cause 
and effect. Both appear to be effects produced by a common cause, namely, 
material ejected from the sun, but different properties of this material seem 
to lead to the magnetic and cosmic ray variations. 

From the study of the 22 Forbush-type events listed in Table I the following 
conclusions may be drawn concerning the relationship between the onset 
times of the decreases and the storm sudden commencements (S.S.C.) for the 
period under consideration: 

(a) The S.S.C. occurred before the intensity was depressed to the east 
of the sun-earth line. 

(6) On some occasions the cosmic ray intensity was depressed from the west 
several hours before the S.S.C. occurred. On other occasions the S.S.C. pre- 
ceded the cosmic ray effects to the west. 


DISCUSSION 


The results presented above indicate that many, if not all, Forbush-type 
decreases conform to a pattern when we take into account the direction of 
maximum sensitivity of a recorder and the onset time observed by it. This 
result seems to arise because of a spatial asymmetry which develops in the 
decrease mechanism at these times. 

It is not yet clear why the decrease should be observed first in directions 
west of the sun-earth line and later to the east. In this connection it should 
be emphasized that the procedure used to determine the directions of maximum 
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sensitivity of a neutron monitor involves the assumption that the conventional 
geomagnetic field configuration is not seriously in error, and that the agency 
responsible for the decreases lies outside the effective field region. If either of 
these assumptions should prove unjustified the directions of maximum 
sensitivity would be different. 

It has been noted that the idealized time-direction diagram for Forbush- 
type decreases predicts characteristic times for intensity maxima that may 
be observed at different stations during the progress of an event. These 
characteristic times correspond rather closely with the time of maximum of 
the daily variation observed at the stations. For example, the characteristic 
time predicted from Fig. 5 for Hobart is 2300-0400 U.T., which is to be com- 
pared with the time of maximum of the average daily variation at Hobart 
during 1957, namely about 0200 U.T. An enhanced daily variation during 
Forbush-type decreases and geomagnetically disturbed periods has been 
reported by others (Elliot 1952). It is possible that the quiet-day daily varia- 
tion may be due to an asymmetrical mechanism similar to that responsible 
for the intensity fluctuations during Forbush decreases. The continuous flow 
of matter from the sun, or solar ‘‘wind”’ (Van de Hulst 1953) may be able to 
set up an intensity-controlling mechanism similar to that due to matter 
ejected from active solar regions. 

The differences in onset times of Forbush-type decreases observed at 
different stations provides yet another piece of evidence that there is not a 
one-to-one correspondence between the decreases and magnetic storm sudden 
commencements, for the latter are known to coincide to within 1 minute at 
widely separated observatories. The detection of cosmic ray changes on some 
occasions for particles arriving from the west of the earth-sun line several 
hours before the S.S.C. suggests that the agency responsible for the decreases 
is able to influence the intensity at a considerable distance from the earth. 
If this is so it means that at least part of the Forbush decrease mechanism is 
located far from the earth and is not geocentric as suggested by Parker (1956). 

Since the cloud of solar material travels at about 5X 10° km/hour, and since 
some preliminary decreases lead the S.S.C. by about 10 hours, this distance 
can apparently be as great as 50X10° km (0.33 A.U.). However, conclusions 
based on the time relations between cosmic ray decreases and S.S.C.’s must 
be regarded with caution because different features of the ejected solar 
material may be responsible for these two phenomena (this is suggested by 
the occurrence of some S.S.C.’s without associated cosmic ray intensity 
changes and vice versa). It does not seem possible to explain cosmic ray effects 
originating some distance from the earth in terms of deflection by magnetic 
fields frozen into clouds of ejected solar matter—in such fields the original 
isotropy of the radiation is expected to be preserved. If anything, the fact 
that the clouds are approaching the earth will result in the energy of the 
albedo from the cloud being increased, hence an increase in intensity would 
be produced. Electric fields in polarized clouds of ionized matter ejected from 
the sun offer a possible explanation. However, there appear to be so many 
inadequately determined parameters in the study of cosmic ray intensity 
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fluctuations that it seems scarcely worth while at present to attempt a quanti- 
tative analysis to test the various mechanisms that have been proposed in 
recent years. 
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APPROXIMATE WAVE FUNCTIONS OF Pb+++ BY THE 
METHOD OF SELF-CONSISTENT FIELD WITHOUT EXCHANGE! 


J. F. Hart AND BEATRICE H. WoRSLEY 


ABSTRACT 


The FERUT program previously described for calculating Hartree-Fock wave 
functions by the method of the self-consistent field has been adapted to the 
configuration Pb*+*+*. Although the exchange factors were omitted, the pro- 
gram was extended beyond its original scope in other respects, and an assessment 
of the difficulties so encountered is made. It might be noted, however, that, 
except in the case of the 4f wave function, it was possible to begin all the integra- 
tions at a common point. Initial estimates were made from the Douglas, Hartree, 
and Runciman results for thallium. The estimates for the core functions were not 
assumed to be satisfactory. The errors in the final wave functions are considered 
to be no more than one or two units in the second decimal digit. 


1. INTRODUCTION 

The FERUT program for the calculation of atomic structures (Worsley 
1958) was designed primarily to obtain wave functions with exchange for the 
atoms up to the end of the second long period. However, with a slight 
modification, the program could in principle be used to calculate wave functions 
and energies of any atom in the periodic table, provided the exchange terms 
were neglected. In this paper we present the self-consistent field calculations 
for the Pb+++ ion as computed by the modified program, and we outline the 
difficulties encountered. 

The neglect of exchange is dictated by the prohibitive amount of machine 
time which would be required on a medium-speed computer. However, this 
simplification seems reasonable for heavy atoms when the relatively more 
important correlation and relativistic effects are being neglected. 

The interest in Pb+++ arises from a spectroscopic study of the ions of lead 
being made at the McLennan Laboratory under the direction of Professor 
M. F. Crawford. Results of this investigation will be published later, but in 
view of the many uses of atomic wave functions and in view of the value of the 
tables as an aid to those engaged in computer calculation, it was decided to 
publish the wave functions separately. Since the functions use the same 
independent variable as used by Ridley (1957) for uranium, namely, 
p = In 1000r, they will be particularly useful for interpolation. 


2. DISCUSSION 
In the case of Pb+++ the Hartree equations to be solved may be written in 
the form 
(2.1) Pi (r) +8i(r). Pir) = 0, 
where 1 = 1, 2,3... 13, 
1Manuscript received April 30, 1959. 
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subject to the boundary conditions 
(2.2) P,(~) = 0, 


and the normalization conditions 


(2.3) f Pi(r)dr = 1. 
0 


In the exchange calculations there are two energy parameters in each 
differential equation so that both of the above conditions are required, whereas 
in the calculations without exchange the condition at infinity is sufficient to 
determine each wave function to within a multiplicative constant. Leaving 
this constant arbitrary until the end of a calculation permits it to be used as 
a scaling factor to prevent overflow during the numerical integration of 
equation (2.1). 

The choice of a suitable numerical method and interval size becomes 
critically important for heavy atoms. The Numerov method, with an accuracy 
of five decimal figures and a speed of two steps per second, was considered 
satisfactory, although the speed could perhaps be improved. On the other hand, 
the choice of logarithmic independent variable, while it obviates the necessity 
of jumps in the interval of integration, slows the calculation appreciably by 
requiring too many points in the region of the core. For example, in the 
integration of the 6s wave function of Pb+++, more than 50% of the points of 
integration lie within r~1/10 atomic units although the function extends to 
r~7 atomic units. If a number of structures were to be calculated with a 
millisecond computer such as FERUT it might be advisable to start some of 
the integrations further out, or to make use of the fact that slight changes in 
the energy parameter only affect the outer portion of the wave function. To 
simplify the coding, the FERUT program was designed so that, for each wave 
function, the following procedure was used. The starting values were calculated 
from the series given by Hartree (1957, p. 81) atr = 0.001 and r = 0.001105, 
and integration was continued outward past the last node until a test could be 
made of the condition (2.2). The use of a common starting point offers no 
difficulty if the computer has fast floating-point arithmetic. However, since 
the starting values are of the order of r't! there may be a complete loss of 
significant figures when using fixed-point arithmetic. In practice it was found 
that, if full use was made of the 40 bit word length, all the wave functions of 
Pb*** could be scaled satisfactorily except the 4f function, which was con- 
sequently integrated separately using a slower floating-point routine. This 
work on a heavy atom confirms our opinion that, for a millisecond computer, 
the automatization must stop short of allowing the entire calculation to con- 
tinue in the machine without operator interference. It was found convenient 
to improve a single wave function during one production run. The choice ot 
wave function to be improved and the trial energy value was left to the 


discretion of the operator. kor microsecond computing machines it is un- 
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doubtedly possible to code a criterion of self-consistency so that production 
will continue until all wave functions are obtained to a predetermined 
accuracy. 


3. RESULTS 


The initial estimates of the wave functions of Pb*+*+* were obtained from 
the results of Douglas, Hartree, and Runciman (1955) for Au and Tl*. The 
thallium wave functions had been obtained by “pure scaling’’ from Au 
(Hartree 1957, p. 116) for all but the 5d and 6s groups which were integrated. 
In the case of Pb+*+*+ equation (2.1) was integrated for the wave functions as 
a check on the estimates, and it was found that improvement was required for 
a number of inner groups. The procedure adopted was to work outward from 
the ls group to the 5p group, improving only those functions for which the 
values differed by more than 2% from the estimated values. Where improve- 
ment was required, the final wave function was substituted for the estimated 
wave function. In this way it was found that the 3d, 4d, 4f, 5s, and 5p groups had 
to be improved, although only the 5p group was sufficiently in error to produce 
a serious change in the total effective field. Finally the 5d and 6s groups were 
integrated until self-consistent. The entire calculation required 40 hours’ 
production on FERUT. The penultimate and ultimate wave functions differ 
by no more than one or two units in the second decimal, or by about 3°%. The 
“‘ultimate’’ wave functions are given in Table I (pp. 986-988). 

ACKNOWLEDGMENTS 
We are grateful to Professors Watson and Gotlieb for making available the 


facilities of the Computation Centre, and also to the late Professor D. R. 
Hartree for his suggestions and encouragement in the earlier part of the work. 


REFERENCES 


Dovaevas, A. S., Hartree, D. R., and Runciman, W. A. 1955. Proc. Cambridge Phil. Soc 
51, 486. 

HartTrEE, D. R. 1957. The calculation of atomic structures (John Wiley & Sons, Inc., 
New York). 

Riptey, E. C. 1957. Proc. Roy. Soc. A, 243, 422. 

Wors.ey, B. H. 1958. Can J. Phys. 36, 289. 








60 '0 20'0 8t'0 110 00°0 S00 oF 0 £60 O10 98 0 9FO LL'I itt FO 9 FLISLO 0 
s0 0 10°0 910 bl 0 00°0 600 2¢°0 OF 0 80°0 92°0 6¢ 0 c¢c'l 98 T 66 9 SFF9OLO 0 
90°0 10°0 bl O 9t 0 00°0 £0°0 ee 0 ct "0O 90° 0 19°0 12°0 cel O9°T is'9 O88F 100 
90°0 100 Zl 0 810 00°0 20 °0 6% 0 : c0 0 8¢ 0 92°0 LT I ZT i¢'9 P9FELO O 


























40°0 10°0 110 61'0 00°0 z0'0 GZ 0 oF 0 ¢0°0 0S 0 180 10°1 LT zo 9 Z8IZ10°0 
10°0 00°0 60°0 020 00°0 10°0 120 tt 0 €0°0 ¢r'O 980 28°0 161 ze'9 €Z0110'0 
40°0 00°0 80°0 1Z'0 00°0 10°0 810 CFO z0'0 1¢°0 060 “L'0 661 ze 9 #166000 
140°0 00°0 40°0 1Z°0 00°0 10°0 910 cr 0 z0'0 Ig 0 16°0 £9°0 861 L1¢°9 6206000 
40°0 00°0 90'°0 1Z'0 00°0 10°0 €1'0 cr 0 10°0 120 06°0 ¢¢°0 161 819 991800 °0 
ms 40°0 00°0 c0'0 020 00°0 00°0 It0 tr 0 10°0 £2'0 880 cr 0 18°1 c6'S 68¢200°0 
o 
= 10°0 00°0 +0'0 020 0c’ 0 00°0 60°0 €F'0 10°0 61°0 cso 80 tert 02'S 9899000 
ns 90'0 00°0 £0°0 610 00°0 00°0 800 It 0 10°0 910 z8 0 oe 0 GL Gr'S 009000 
: 90'0 00°0 €0°0 810 00°0 00°0 10°0 6£°0 10°90 €1°0 82'0 12'0 L9°1 91°¢ FLFS00'0 
a 90°0 00°0 Z20'0 LV-6 00°0 00°0 90°0 1¢°0 100 110 ¢2°0 Ze 0 6G 1 88 F €66F00°0 
> 90'0 00°0 z0'0 910 00°0 00°0 c0'0 00°0 60°0 0L°0 610 6F I 1o°F Z8FEOO 0 
v3 
2 c0'0 00°0 Z0'0 cl'0 00°0 00°0 +0'0 €£°0 00°0 80°0 19°0 91°0 Ih TI 68°F ¢coron 0 
g ¢0'0 00°0 10'0 #10 00°0 00°0 £00 ig’ 0 00°0 90°0 290 €1°0 Ie 1 66° ¢ 699£00 0 
= ¢0'0 00°0 10°0 10 00°0 00°0 €0°0 620 00°0 c0'0 8¢'0 Il'O Zo I [ee 0ZE£00 0 
i: +0'0 00'°0 100 ZI 0 00°0 00°0 Z0'0 ‘0 00°0 +0'0 #o°0 60°0 cl Ccr'¢é £00800 0 
5 £0'0 00°0 10°0 Z1'0 00°0 00°0 z0'0 0 00°0 +00 0¢ 0 10°0 901 61'¢ 8122000 
~~ 
< #0'0 00°0 10°0 110 00°0 00°0 z0'0 £2°0 00°0 £0°0 9F 0 90'0 16°0 C6 °% 09FZ00 0 
2 €0°0 00°0 10'0 O10 00°0 00°0 10°0 1Z'0 00'°0 €0°0 €r 0 160 oe 922200 0 
2 £0°0 00°0 00°0 60°0 00°0 00°0 10°0 020 00°0 z0 0 6£°0 £0°0 £80 C'S £10Z00'0 
> €0°0 00°0 00°0 80'0 00°0 00°0 10°0 810 00°0 z0'0 9¢°0 £00 110 ie'Z Z28100'0 
€0°0 00°0 00°0 80°0 00°0 00°0 10°0 110 00°0 10°0 ¢£°0 €0°0 02°0 fice 6F9100'0 
<= 
= z0'0 00'°0 00°0 10°0 00°0 00°0 10°0 c1l'0 00°0 10°0 0¢'0 z0°0 90 €6°1 Z6F 1000 
Z z0'0 00°0 00°0 90°0 00°0 00°0 00°0 FIO 00°0 10°0 820 z0'0 09°0 Toa O0SE100°O 
“ Z0'0 00°0 00°0 90°0 00°0 00°0 00°0 £10 00°0 10°0 cZ 0 z0'0 cc '0) 291 122100 0 
- Z20'0 00°0 00°0 c0'0 00°0 00°0 00°0 Z1'0 00°0 00°0 £20 10'0 0¢ 0 6F I GOTTO0 0 
z0'0 00°0 00°0 c0'0 00°0 00°0 00°0 Il'0 00°0 00°0 Ze 0 10°0 cr '0 LE 1 0001000 
(S9)d (PS)d (4S)d (SS)d (fF)d (Phd (Gd (Std (P&)d 4S)d (S8)d (F@)d (Sd (SI)d 4 
StS c9o'¢ LIL 086 6 ZI FE @ St 82S 981 SIZ 9€% 666 Z86 6209 a 
sg PS dg SG If Pr dt Sh PE dg SE IZ SZ S| ju 
"a ve ‘0 = dg ‘40001 uy = d ajqeiaea juapuadapu yun payrjnaye> 
= +44 4d JO UOIZEANSYUOD sQo (PG) AY} 1OJ VBuULYOXA YNOYIIM SUOTOUN] BALM [eIPes pazi];eULION 


oO 


I AHTEaVL 











b €20- L¢:0- ¢90 #£69'0-— Z2g'I 99°I- LT tLI- 8F0O 8¢0- 90 #42000 F0S9E 0 
4 €2:0- 6F0- +9°0 69'0- Zl tS'I— Zit zi- 690 &s0-— &6'0 100 00°0 ogoge 0 
IZ 0- 9F0-— 99:0 69'0-— 6F'I PE I- 2'I L4UI- ¢60 8I'I-— Og't z0'0 10'0- 18862 0 
SI'0- 68°0- &F'0 €b0- ail 90 I- €8'°0 *+¢2°0— O81 €o'I- g9o'l c0'0 £0'0- €FOLZ 0 
s0'0-— Of 0- G20 wz o- Zi GLO- OF 0 920- 6F'I ss I- 02 800 #£460°0- 69FFZ ‘0 
00:0 610- 900 000 Oat If'0- 200- 610 8g2Z'I 61 Z—- See 10 8s1'0- 9041220 
400 200- Il'0- 02:0 40°1 10°;0- If 0- 69°0 #02 car2—- SSZ Go  Zz'0- L€€002 0 
Z1°0 100 90- 98:0 60 0 jZ0- 880 992 982- 19% 680 2:0- ZLZI8T'O 
SIO 00 28°0- &SF'0 180 OSO0 6'°0-— +#0'I €h2 982-— 08'% 090 #62'°0- ZZ0F9I'O 
: 910 170 8 &F0- 8F0 89°0 120 g0I- L0°T gS'Z th Z— #22 £&$8°0 fl I E1F8FT 0 
+ 
2 9L'0 170 8 &bF0- 9F0 980 980 0I- 660 98% O062- ZI 911 €o°1 00°0 O6ZFEL'O 
te €1'0 1¢°0 Ir 0-— 880 9F0 G60 G&0- 180 ts GsI- gel Zot 96 1- 20°0 OISIZI°0 
° O10 2:0 Ss0- 2420 2eo- oo! 6L:0- 249°0 &%S &FI- &2°0 161 ZeZ— +0°0 L€6601 0 
Y c00 €80 9%0- SIO 620 66:0 680- O£0 822 $%L°0- S10 Ie'% 148°%2-— 20°0 F8F660 0 
S 10°O0 20 210- 20 0 9%%0 j%Zeo- ¢0'0 Of 6F0-— zeo- 0212 GZ ¢g—- O1'0 110060°0 
me 
> €0'0-— 0f0 =2400-— 60'°0— 810 680 #Ft1'0-— 2 0- O61 OF 0- OL0- SOE Ge'E- GIO ISFI80'0 
5 900- 820 20 #42361'0- F1'0 30 200 #£&F0- 69 1¢°0 IOI- 98'¢ €2¢- O£°0 00L¢20'0 
Fs 60'0-— &2'0 110 92°0- O10 20 92:0 #£69'0-— 8 €2:0 S82 I- 6¢9°¢ 6L:¢- 9F'0 989990 0 
Q Il'O0- 20 £810 I¢0- 800 +90 20 £69'0-— 221 zO'T GhiI- G2¢ 9L¢- 219°0 0F£090 0 
< Il'O0- 61:0 & 0 g¢80- 900 g¢0 #&¢°0 2'0- 80°T #21 SFI- Es ¢ #9°E-— O0'T S6SFS0 0 
> Il0- 91:0 220 &&0- 00 oOF0 #90 &€20- 1160 ~~ 6&'I L4eI- 8 9@E- SEI ZOFGEO 0 
a O1'0-— &s1'0 62:0 I 0- £00 680 690 89':0- 92:0 &8F'T €eI- 22°¢ 48°2- =I8'T 1OLFFO'0 
B O1'0- I1'0 I¢'0 s20- 200 zs0 22°0 190- 29°0 eI 02 I- 98 @bZ- #2 LFF0FO 0 
3 80'0- 60°0 ig¢0 €0- 200 90 #£22°0 if 0O- 190 OG'T 660- LE 1461- 822% 866920 0 
2 900- 2400 Of0 #£8t0- 10:0 70° 86©02:0)»=668'0- =O«IF'0sCO*FNF'T 9L0- 92'¢ 6EI- ees CLTEso 0 
s 
Z 410 0- 900 62:0 ¢10- 100 200 #290 %L4%0- £80 ~~ 8E't €¢0- ZOE 16'0—- 98'°¢ £96620 0 
be 20'0- 00 242°0 iS o=- Wo FLO £90 Si g=— v6 0 62 I 660—- 926 OF 0- OF F E1TZZ0°0 
x 10 0- 00 G20 £2 0- 000 i100 }§=6 880 t00- 120 61 100- 18% 900 68°F ELSFZO'O 
= 100 00 €%0 #%$£00 oOO 600 €¢0 £900 910 ~~ 801 €l'0 ZZ ce ae ee SO 8612200 
200 8=—00 120 10 OOO 200 gO Sto €t0 £2460 20 10°% Ig0 »=6 92 9800200 
(s9)¢d (PS)d (4S)d (88)d (hd (Prd (4d (hd (Ped Wd (SE)d Bd CB@d  (SDd 4 
srz.6S9 COLT LO 6COCG'T +1 oSbh $s 29 981 81Z 9% GEG Z86 6209 My 
sg PS dg s¢ ff Pr qt Sh PE d¢; S§ qdZ SZ ST ju 


LO = 92 ‘4Q001 Uy = 9 a]qeiaea Quapuadapul YIM pazEpNoyeD 
444d JO UONRINSYUOD S$Oo(PG) BY JO} 9BuULYOX NOY SUOTIUN, 9AVM [LIPe4 pazIPeULION 


(panunuod) | ATAVL 









CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


988 


ge 
20° 0- 
40 0- 
80°0- 
cE U= 


02° 0- 
62 0- 
oF 0- 
1$°0- 
19°0— 


69°0—- 
yL 0-— 
pL 0—- 
. 0- 
€9°0- 


og 0- 
L¢°0- 
12'0- 
go: 0- 
oro 


£2'0 
oe 0 
9¢°0 
9¢°0 
0e 0 


co 0 
Il'0 
1 o= 
EE-@= 
61° 0- 


(S9)d 


8h 3 
Sg 


So Q “4 
a=S35 558 


ESSERE CITAA = , 
coooo Sooeoo coooo ofFeocoo SoS 


Crwws9s AN AnowrnN 


mud tH 
Hin OO 


g9°€ 
Ps 





*31Qe} SIY? UL pasn are SPUN S1WOZe VaITIeH{ ALON 





L61E8 2 
FSFE’ 
16Z00°9 
991gh'S 
LL¥I6'F 
LOLEE FE 
18€20'# 

00°0 S609 '€ 
10°0-— 00°0 LEP6Z E 
€0':0- 10°0 96086 Z 
c0'0- 20°0 82269 'Z 
60°0- ¢0'0 00°0 O90FF Z 
ST'0- 60°0 10°0 SE80Z Z 
t20- SIO 200 02866 I 
ce'0-— &¢2'0 ¢0'°0 #00°0 $0808 T 
6F'0-— 28:0 #¢0°0 10;0-— 00°0 86629 'T 
79°0-— €9°0 800 z200-— 100 #£00°0 OL08F I 
os'0- 120 &10 #£=+0'0- 20:0 10°0- Er6eE I 
96'0- 680 810 g0°0- #00 #££€0°0- 16113 'T 
s0°'I- 90°T cZ'0 36 ST'0-— +=80'0 = 90 0— €9960'T 
CI'I- 8st GeO 1Z0- SIO 21'0- 122660 
9I'I- eT 4vO 806 88'0- ~«930 2 = 0- 68268 0 
80'I- OZ T 09°0 «=68FO0- OF0 98°0- IFZI8 ‘0 
Z6°0- 90'°T c¢40 249°0- 090 #£=98'0- 00:0 000 #4000 OISE2 ‘0 
0L0- #80 O60 £88°0- +80 1g8'0- 10°0 100- 10°0 FIS99 0 
€F0- 9°00 SOT II- Ot s0I- ¢€0°0 ¢0°0- ¢€0°0 SgT09 ‘0 
FL O0- 2:0 ~~ 61'T ceiI-— sel 98 I- 900 90°0- 90°0 LStFS ‘0 
FLO Il'O0- IT IG'I— sgt 09 I- If'O0 2Zt'o- &1'0 C1Z6F 0 
se'0  s80- Zl 9 I- ELT LLI- 61'0 %'0-— +20 O8SFF'O 
cg¢'0 Ss 8S O—_~=sCéGF'T 69'I- 82'I csI- 1€°0 148°0- I¥'0 €FE0r 0 
(4q)q (sc)d (Std (Phd (4d (Std (Ped (F8)d (Sed (Fd (Std (Sd 4 

LE 3 086 °6S1 rIe 2h 82g 98T 81Z 9% 646 286 6209 1M9 

dg sg St PY qt SH PE dg Sg qZ SZ ST 1u 








LO = 99 ‘40001 Ul = 9 a[qeizea Juapuadapul 4M parejnaje> 
444d JO UOIZRINSYUOD $Qo,(PG) 9Y} 10} aBULYOXa JNOYIM SUOI}IUNJ BAM [LIPeI Pazi[euION] 


(papnjau0D) | ATAVL 








A CRYOSTAT FOR MEASURING SPECIFIC HEATS BETWEEN 
1 AND 4.2° K USING A MECHANICAL THERMAL SWITCH! 


F. D. MANCHESTER? 


ABSTRACT 


A design for a low temperature cryostat using a mechanical thermal switch 
is described and its performance and some general features of interest in its 
operation, such as the importance of vibration, are discussed. Sample measure- 
ments on pure copper and a cupro-nickel alloy are presented. 


In recent years interest has been revived in a design of cryostat for specific 
heat measurements which makes use of a mechanical thermal switch instead 
of helium exchange gas for cooling the specimen to the working temperature. 
This method, first tried out by Keesom and Kok (1933), has been used in 
cryostats designed by Ruehrwein et al. (1943), Westrum et al. (1953), 
Kostryokova and Strelkov (1953), Dolecek and Madden (1953), Ramanathan 
and Srinivasan (1955), Webb and Wilks (1955), and Rayne (1956). A discussion 
of the undesirable effects produced by exchange gas has been given by Webb 
and Wilks (1955) and Corak et a/. (1955) and their occurrence in experiments 
has been described by Keesom and Kok (1933) and Horowitz et al. (1952). The 
present paper describes a cryostat using a mechanical thermal switch patterned 
after the design of Webb and Wilks and presents some information and 
experience gained from its use in specific heat measurements. This design was 
chosen because it seemed to offer the best scope for extension of its application 
to other experimental situations, after some experience had been gained in 
using it for calorimetry. 


DESCRIPTION OF THE APPARATUS 


The general layout of the mechanical thermal switch cryostat is shown in 
Fig. 1 and the details of the thermal switch in Fig. 2. 

The calorimeter C was suspended by nylon threads from the frame F and 
the space surrounding the calorimeter was maintained at a pressure of about 
10-§ mm Hg by continuous pumping through one of the lines P. Temperatures 
between 1.1 and 4.2° K were obtained by pumping on the liquid in the inner 
liquid helium chamber H and the thermal switch was used to make contact 
between the chamber bottom and the sample. The mechanism M of the thermal 
switch was controlled by the actuating rod R, which passed through a tube to 
the room temperature end of the apparatus, where it left the vacuum system 
through a flexible bellows S and connected with a screw arrangement used for 
raising or lowering the rod. The switching action was such that the rod was 
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Fic. 1. General arrangement of mechanical thermal switch cryostat. 


under tension when the switch was closing, i.e. when the maximum force was 
required. A total force of up to 6 kilograms could be applied between the jaws 
of the switch. The rod was thermally anchored at three points A, along its 
length. One, just under the Dewar cap D, was just below the temperature of 
liquid nitrogen; another, just above the main part of the cryostat, was at a 
temperature of 4.2° K; and the third was at the bottom of the rod at a tempera- 
ture corresponding to that of the pumped inner helium chamber. The switch 
mechanism consisted of the copper jaws J, faced with a layer of indium metal 
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Fic. 2. Detail of mechanical thermal switch. 


approximately 1 mm thick, which could be made to grip the cylindrical 
clamping post K. Motion was imparted to the jaws through the switch arms 
moving about the pivots V and coupled to the Teflon head on the end of the 
actuating rod. A combined extra radiation shield and guide for the actuating 
rod was provided at G and a copper disk X was used as a further radiation 
shield and a thermal anchoring point. Teflon was used for the head and for 
bushes in the pivots V, to minimize any heat produced by friction between the 
moving parts of the switch. Copper braid B was used to connect the jaws of 
the thermal switch to the copper posts N on the bottom of the inner helium 
chamber. This provided a direct connection which did not depend on the com- 
ponents of the switch mechanism. The clamping post K was a small copper 
cylinder drilled out to reduce its mass and was well polished on its outside 
surface. Initially the jaws which engaged with the post were also of polished 
copper but were later faced with indium metal. The difference in performance 
which resulted can be seen in Fig. 3. These indium-coated jaws proved quite 
serviceable in use, no trouble was experienced with the indium sticking to the 
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Fic. 3. Comparison of cooling rates obtained with helium exchange gas, indium-clad switch 
jaws, and copper jaws. 


clamping post, and the jaws remained usable with no special treatment, other 
than working the switch a few times to improve the matching of the surfaces, 
after repositioning the calorimeter. Following Webb and Wilks an insulated 
wire was connected to the calorimeter and a simple indicator circuit arranged 
to tell when the thermal switch had made contact with the calorimeter. The 
electrical leads to the calorimeter were brought through the tubing shown in 
Fig. 1 to the seal S, then into the helium bath at 4.2° K, and from there to 
room temperature. 


OPERATION OF THE MECHANICAL THERMAL SWITCH 


The calorimeter was suspended in the frame by its nylon threads and care- 
fully adjusted to be centrally placed with respect to the jaws of the switch, so 
that no discernible movement of the calorimeter occurred when the switch 
was worked. 

Cooling with the switch was normally done in three stages: room tempera- 
ture to liquid nitrogen temperatures, then to 4.2° K, and then down to as low 
as 1.1° K. In many cases the specific heat of the sample was large enough that 
the switch had to be connected only lightly when cooling from liquid nitrogen 
to liquid helium temperatures to prevent excessive boil-off of liquid helium. 
In cases where this was not a limitation (for example, for 75 g of copper) a 
sample could be cooled from room temperature to below 2° K in less than 4 
hours. The force applied to the switch was also determined by the temperature 
at which it was intended to be opened and the thermal capacity of the body 
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being cooled. The larger the force applied, the larger the heating produced at 
the parting surfaces when the switch was disconnected, and the resultant 
temperature rise of the cooled sample would depend on its thermal capacity. 
At 4° K the specific heat of the sample was in most cases sufficiently high that 
the amount of heat developed at the break was not important. Thus down to 
4° K the force on the jaws could be large, and particularly from room 
temperature to liquid nitrogen temperatures it could be a maximum. In cool- 
ing from 4° K down, however, a balance had to be struck between a loading 
sufficient to cool the sample at a reasonable rate and one not so great that the 
heat developed at the break was inconveniently large. This could mean that 
one had to wait up to 1 hour to cool down to about 1.1° K but if one did this, 
the switch could be disconnected with almost no warming up at all (see Table I). 





TABLE I 
Thermal data on performance of mechanical switch 
Sample run No. i, i Ars Q, joule 
1 1.124 0.004 1.221075 
2 1.163 0.002 0.851075 
2 1.145 0.005 2.09X10-5 
3 1.143 0.001 0.64X10-5 
6 1.187 0.002 1.29X1075 
Webb and Wilks (1955) — = 10-¢ 
Dolecek and Madden (1953) — — 0.75-2.25x10-4 


Apart from the force applied to the jaws of the switch the other most 
important consideration was vibration, that is, the mechanical vibration 
transmitted through the framework of the apparatus to the switch mechanism 
and the sample. The role of vibration as a source of extraneous heating in low 
temperature experiments employing thermally isolated samples, such as in 
calorimetry and adiabatic demagnetization work, is already widely appreciated 
(Hull, Wilkinson, and Wilks 1951; Darby e¢ al. 1951; Wheatley, Griffing, and 
Estle 1956). However, in the case of the mechanical thermal switch the vibra- 
tion also affects the action of the switch itself. Firstly, the amount of extraneous 
heat received by the sample and the thermal conductance of the switch deter- 
mine how nearly the sample can be cooled to the temperature of the pumped 
helium bath, and vibration is a very important source of extraneous heating. 
Secondly, the amount of heat produced at the disconnecting of the switch is 
affected, because the relative vibration of the jaws and the clamping post 
produce a heating additional to that which results just from the friction of the 
surfaces parting simply. For these reasons it is essential to have vibration in the 
apparatus reduced as far as possible; indeed, the extent to which this can be 
done essentially determines the useful limiting temperature to which the switch 
can cool the sample. In the case of the present apparatus the site on which it 
was placed initially was such that the extraneous heating of the sample pro- 
duced by vibration was more than 200 ergs/minute. Similar figures have been 
quoted by Hull, Wilkinson, and Wilks (1951) and Darby et al. (1951). Vibration 
of this magnitude made calorimetric measurements with metals virtually im- 
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possible and the performance of the thermal switch quite misleading. Collins 
and Zimmerman (1953) refer to the importance of heat generated by vibration 
in their description of a mechanical thermal switch used in a two-stage adiabatic 
demagnetization apparatus, although not specifically to its effect on the switch- 
ing action. Once the apparatus had been placed on a relatively vibration-free 
support the extraneous heat input to the sample was reduced to a tolerable level 
and the performance of the switch was much improved. Even so, vibration in 
the apparatus varied from day to day and with the time of day, so that in 
general, except for quiet days in week ends, it was not possible to operate the 
thermal switch, during the day, to cool a sample to the lowest parts of the 1 to 
4° K temperature region and disconnect it satisfactorily. All work in this region 
had to be left until late evening. 


THE PERFORMANCE OF THE MECHANICAL THERMAL SWITCH 
In Fig. 3 comparative cooling rate curves are given for a sample of 75 g of 
copper being cooled from about liquid nitrogen temperatures to the boiling 
point of helium. It can be seen from this that the cooling rate with the thermal 
switch is sufficiently close to that obtained with exchange gas that there is no 
loss in experimental convenience. In Fig. 4 the temperature dependence of the 
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Fic. 4. Temperature dependence of the thermal conductivity of the switch in the region 
below liquid nitrogen temperatures. 


thermal conductance of the thermal switch (using indium-clad jaws) is shown 
for the region between liquid nitrogen and liquid helium temperatures and in 
Fig. 5 it is shown for temperatures between 4 and 1° K. The temperature 
dependence of the conductance is similar to that found by Berman (1956), 
close to 7? in the 2-4° K region and gradually decreasing to become essentially 
temperature independent at liquid nitrogen temperatures. These results are 
not presented as a continuous curve as they were taken while routine cooling 
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Fic. 5. Temperature dependence of the thermal conductivity of the switch in the liquid 
helium temperature region, 


was being performed, with different forces on the jaws in the two temperature 
regions. Nevertheless, the curves in Figs. 4 and 5 represent the typical tem- 
perature dependence and magnitude of the thermal conductance. The figure 
quoted by Webb and Wilks (1955) for the thermal conductance of their switch 
is ‘of the order of ’’ 1.5X 10-4 joule/deg sec, presumably in the liquid helium 
temperature range. Dolecek and Madden (1953) found a temperature depen- 
dence of 7'*X 10~° joule/sec deg for their switch in the liquid helium tempera- . 
ture range. This gives a value of 7X10-° joule/sec deg at 4° K, which is 30% 
less than the value obtained from Fig. 4. Greater conductances could be 
obtained in the liquid helium region by applying more force to the switch but 
the values given in Fig. 5 are those obtained from experiments in which the 
switch was successfully disconnected with very little thermal disturbance. 
This would not be the case if the switch was clamped tightly as was discussed 
earlier. 

Table I gives typical data on the warming up, Q, produced in the sample by 
the opening of the switch. The data are taken from observations made when 
vibration in the sample was very small and represent the warming produced 
essentially by the parting surfaces. No discernible warming was produced when 
the switch was closed, which is in contrast to the findings of Dolecek and 
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Madden (1953). When vibration conditions were bad the warming would be a 
good deal greater and the temperature drift of the specimen after disconnecting 
the switch made it difficult to tell how much warming had been produced just 
by the switching action. The table also shows typical temperatures to which 
the sample was cooled and the useful temperature at which measurements 
could be started, i.e. T+AT7. That measurements can be started at such 
temperatures is one advantage of a thermal switch; as in cases where the heat 
input due to vibration is reasonable, specific heat measurements can be started 
immediately the switch is disconnected. When exchange gas is used there is a 
delay while the exchange gas is pumped out to a sufficient extent for measure- 
ments to begin and in this period the sample will generally warm up, as the 
temperature drift is usually in the direction of warming at the lowest tempera- 
tures. Also, if it is desired to repeat measurements at the low temperature end, 
this can be done very conveniently using the switch. 

The thermal isolation of the specimen was very good and for most of the 
temperature range, the drift of the sample temperature was compatible with 
the conduction of heat, along the nylon suspension threads, to the helium being 
pumped in the chamber H. It was usually only necessary to raise the tempera- 
ture of the pumped helium once in the course of taking measurements between 
1 and 4° K and for some substances with large specific heats there was no need 
at all. If vibration had been less of a problem, then finer suspension threads 
could have been used and the heat leak further reduced. 

The foregoing represents the essential data on the performance of the 
mechanical thermal switch. These observations were made in the course of 
using the apparatus for measurements of specific heat and were to some extent 
incidental to the main purpose of the experiment. Thus no measurements were 
made of the force applied to the jaws during the cooling down and of the 
corresponding thermal conductivity. Some low temperature measurements of 
thermal conduction across metallic surfaces in contact and under varying loads 
have been made by Jacobs and Starr (1939) and more recently by Berman 
(1956). These studies are of help in understanding the processes involved and 
in designing a mechanical contact, although more work remains to be done 
in this field. 


SPECIFIC HEAT MEASUREMENT 


Experimental Details 

Temperatures were measured with a 10-ohm, 1/2-watt, Allen Bradley 
resistor calibrated against the vapor pressure of liquid helium. The two- 
constant interpolation formula (Clement et al. 1955a) 


V (log R/T) = a log R+b 


was used in the manner suggested by Keesom and Pearlman (1956). The 
thermometer was calibrated by admitting exchange gas to the calorimeter 
space at the conclusion of a set of specific heat measurements on a sample and 
a total of about 24 calibration points taken over the range 1.1 to 4.2° K. This 
calibration procedure had to be repeated for each run made at liquid helium 
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temperatures, after the cryostat had been allowed to return to room tempera- 
ture, as under these circumstances the calibration curve was displaced in an 
unpredictable fashion from run to run. As an example of the magnitude of this 
displacement, the two most extreme cases out of seven runs with one particular 
thermometer showed a change of 0.050° K for the temperature corresponding 
to a resistance of 100 ohms. Vapor pressure readings were corrected where 
necessary for thermomolecular pressure difference using the tables of Roberts 
and Sydoriak (1956) and converted to absolute temperatures using the 1955 E 
scale (Clement et al. 1955b). The resistance of the carbon resistor was measured 
using a method developed by Dauphinee and described by Dauphinee and 
Mooser (1955). This method has the advantage that thermoelectric e.m.f.’s 
had no effect on the measurement and that the thermometer resistance was 
compared directly with a reference resistance decade and the comparison 
presented continuously on a chart recorder. The temperature drift before and 
after the heating period was measured from the chart and then extrapolated 
to the middle of the heating period, as described by Keesom and Kok (1932). 
The temperature drift varied from zero to 10-°° K/sec. Thermometer currents 
varying from 2 wa to 15 wa were used over the temperature range, and within 
the required precision of the thermometry, this variation of the current had no 
effect on the measured resistance. The size of the temperature increment used, 
AT, was kept such that the resultant error in the measured specific heat 
{4(AT/T)? for C « T*} was less than 0.1% for all parts of the temperature 
range. Errors arising from the calibration of the thermometers and from the 
vapor pressure — temperature relations themselves have been discussed by 
Corak et al. (1955) and by Rayne (1956). The calibration procedure in the 
present work has been carried out in such a way that the errors involved are 
well within the limits discussed by these authors. On some occasions when a 
specific heat measurement was being made, there were small amounts of energy 
fed into the heater circuit from induction effects produced by the shutting off 
of the neighboring electric clock circuit. These stray energy inputs were largely 
eliminated but could be noticed now and again near 1° K. As, however, they 
showed up plainly on the chart record, suitable correction could easily be made. 

The calorimeter heater was wound from manganin wire and had a resistance 
of about 2000 ohms. The electrical leads connecting the heater on the calori- 
meter with the rest of the circuit on the body of the cryostat were made from 
manganin wire with an outer coating of soft solder. The soft solder became 
superconducting at liquid helium temperatures, providing leads which com- 
bined zero electrical resistance with high thermal resistance. The power 
dissipated in the heater was measured using a potentiometer to measure the 
current through, and the potential drop across, the heater, in conjunction 
with an electric clock, driven by a standard 60-cycle supply, which was used 
to measure the duration of the heating period. 

The calorimeter arrangement depended on the substance being measured. 
For most of the measurements, which were on copper alloys, a “‘jacket’”’ type 
of calorimeter was used. This consisted of a copper cylinder 1} in. long provided 
with a slightly off-center hole of 1.000-in. diameter along its length. Into this 
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hole the carefully machined sample could be slid and clamping screws used to 
tighten the jacket around the sample. Clamping post, thermometer, and 
heater were permanently attached to the jacket and thus once the specific heat 
of the jacket had been measured, successive samples could be measured with a 
minimum of trouble. The heater was wound on a small copper bobbin placed 
in a hole running through the length of the jacket wall and the bobbin was 
soldered to the jacket at one end of the hole. The carbon resistor had its outer 
plastic covering ground off and after being covered with an insulating layer of 
baking varnish it was slipped into another hole in the jacket wall and held in 
place with more varnish. The thermal contact between the jacket and the 
sample was good and for all of the samples measured, there was never any 
indication of a lack of thermal equilibrium in the calorimeter. For the measure- 
ment of the specific heat of rubidium (Manchester 1959) the calorimeter was a 
cylindrical copper shell provided with similar fittings. 


Results 
As a check on the performance of the apparatus the specific heat of pure 


copper was measured. The copper, of nominal purity 99.999%, obtained from 
the American Smelting and Refining Company, was melted under vacuum 
using an induction furnace and a graphite crucible. For this sample the values 


y = 0.696+0.021 millijoules/mole deg? 

bo = 344+7° K 
were obtained by using the method of least squares to fit a straight line to the 
data plotted in the form of C/T vs. T*. The error in y and 49, assumed to be 
random, was estimated from this fit at the 95% confidence level. The scatter of 
the measurements is greater in this particular case than in later experiments, 
because this was the first experience with this apparatus. However, the mean 
values for y and 6) compare reasonably well with those quoted by Corak et al. 
(1955), Rayne (1956), Griffel et a/. (1957), as shown in Table II. It is possible 











TABLE II 
Comparison of values for specitic heat of copper 
Oo, °K y millijoules/mole deg? References 
344.543 0.687+0.012 Rayne (1958) 
342+2 0.691+0.006 Griffel et al. (1957) 


34447 0.696+0.021 This work 


that the fact that the copper was melted in a graphite crucible enabled some of 
the metallic impurities to enter into solid solution and thus contribute to a 
higher electronic specific heat. Following this check with copper, measure- 
ments were made of the specific heat of some dilute iron-in-copper alloys. 
These will be discussed in a separate publication. Also a measurement was made 
on a cupro-nickel alloy. This last measurement was in the nature of an explora- 
tory one, undertaken for the following reason. Some measurements on cupro- 
nickel alloys had been made by Keesom and Kurrelmeyer (1940) and from 
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these they obtained a curve giving the dependence of the electronic specific 
heat on the composition of the alloy, a curve that has not been very amenable 
to theoretical interpretation. Because of this and because of the behavior of 
some of the C/T vs. T* plots shown in Keesom and Kurrelmeyer’s paper (a 
pronounced upward displacement of the curve at the low temperature end, see 
also Keesom and Pearlman 1956), it was decided to re-examine the specific 
heat of a cupro-nickel alloy. It seemed most likely that the unusual behavior 
of the C/T vs. T* curves could be due to the helium exchange gas which Keesom 
and Kurrelmeyer had to use, in order to reduce the thermal equilibrium time 
of their calorimeter to a reasonable value. In their estimate of the contribution 
to the measured specific heat made by the helium gas contained in their 
sample, Keesom and Kurrelmeyer only took into account the specific heat of 
the gas and not the heat of its desorption. Actually, the heat of desorption for 
helium gas was not known until it was measured, later on, by Keesom and 
Schweers (1941). The heat of desorption is considerable, however, and would 
be quite capable of accounting for the observed behavior of the measured 
specific heat. 

To investigate this point and to see if measurements could be made without 
using any exchange gas in the sample, a cupro-nickel alloy (80% nickel) was 
prepared and its specific heat measured. The alloy was prepared from spectro- 
scopically pure nickel and copper, melted by induction and chill-cast under 
vacuum in a split copper mold. After casting, the ingot was annealed for 72 
hours at a temperature of 1000° C to remove coring, and was then machined 
to size. Even with the sample prepared in this way, there were variations in 
composition of as much as 3% between ends, as determined by chemical 
analysis. However, for the exploratory measurement undertaken in the present 
instance this is not a serious matter. In the measurement of the specific heat 
no trouble at all was experienced with lack of thermal equilibrium in the sample 
and the measurements were straightforward. This was presumably because of 
the good thermal contact afforded by the jacket design of calorimeter. Thus 
even though the thermal conductivity of samples of 50 to 60% nickel would be 
still less, the experience with the 30% sample is good indication that a full 
range of alloys could be measured without any real difficulty being encountered 
from lack of thermal equilibrium, especially as the thermal isolation of the 
sample was sufficiently good that some lack of thermal equilibrium could be 
tolerated anyway. 

Figure 6 shows C/T vs. T° plotted for the 30% nickel alloy (curve (6)) and 
compared on the same scale with data for two of the alloys ((a) 40% nickel, 
(c) 20% nickel) measured by Keesom and Kurrelmeyer. The ordinates have 
been displaced in order to show the plots more clearly. It can be seen that for 
the 30% alloy the experimental points lie fairly close to a straight line and do 
not show the systematic deviation from it which was found by Keesom and 
Kurrelmeyer for the 20% and particularly for the 40% alloy. Thus probably, 
although not conclusively, this effect was due to the desorption of exchange gas 
inside the sample. Certainly, the measurements on the 30% alloy indicate that 
more reliable measurements of y and 40 could be obtained with apparatus of the 
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Fic. 6. C/T vs. T? for cupro-nickel alloys: 
(a) Cugo Nigo (Keesom and Kurrelmeyer 1940), 
(6) Cuzo Nizo (this work), 
(c) Cuso Nizo (Keesom and Kurrelmeyer 1940). 


All the data have been plotted on the same scale but the ordinates have been displaced for 
clarity. 


type used in the present work. The measured specific heat of the cupro-nickel 
alloy (31.3 at.% Ni, effective atomic weight = 62.04 g) is represented by 

y = 2.027+0.013 millijoules/mole deg?, 

6) = 375+5° K. 


The errors in y and 6) were estimated from straight-line fits to plots of C/T vs. 
T?, evaluated at the 95% level of confidence. 

In a publication which appeared since the completion of the present work 
(Guthrie et al. 1959) further measurements of the specific heat of some cupro- 
nickel alloys have been reported. For copper concentrations of less than about 
75% these authors find a negative slope for C/T vs. T?, which, however, is not 
as discontinuous as that found by Keesom and Kurrelmeyer, and Guthrie et a/. 
attribute this rise in specific heat to a magnetic contribution. However, some 
of these measurements were carried out with helium exchange gas inside the 
calorimeter, performing a similar function to that used in Keesom and Kurrel- 
meyer’s experiments, so that some of the results could still be misleading, 
because of the contribution from helium gas desorption. If the measurements 
on solid bars were made without any internal exchange gas, as seems likely, 
then there would appear to be a genuine effect associated with the specific heat 
of the alloy; but to make interpretation more certain, measurements should be 
made with an apparatus such as is described in the present paper, where 


exchange gas is not necessary. 
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CONCLUSION 


A mechanical thermal switch cryostat has been developed along the lines of 
the design by Webb and Wilks. Its performance and some features of interest 
in its operation have been described; for instance, vibration is a very important 
consideration in determining the performance of the switch. Reasonable and 
practical cooling rates can be achieved without applying a large force to the 
switch and, in fact, a large force is undesirable if warming at the breaking of 
the switch is to be minimized. The switch has been successfully used for cooling 
metal samples in specific heat measurements in an arrangement which should 
be readily applicable to other experimental situations. 
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ORBITAL FOLLOWING IN TETRAHEDRAL MOLECULES'! 


LESSER BLUM AND NorAuH V. COHAN 


ABSTRACT 


On the basis of Heath and Linnett’s formulas, the ‘“‘orbital following”’ force 
constants of NH3;- and NHj-like molecules are calculated. The anomalous 
results are interpreted as a consequence of orbital deformation during the 
vibratory movement. 


For the study of molecular phenomena it is of interest to know what happens 
to the electronic charge distribution during vibration. According to the Born- 
Oppenheimer (1927) approximation the charge is disposed in such a way as 
to minimize the electronic energy. On the basis of similar considerations Heath 
and Linnett (1948a, 6) introduced the concept of “‘orbital valence force field’’: 
they supposed that the function can be taken as resulting from the contribu- 
tions of (a) the interaction of the nuclei with the electronic clouds and (6) the 
interaction between the electronic clouds; and, as the electronic clouds are 
disposed so as to minimize the electronic energy for each nuclear configuration, 
the electronic clouds do not follow this movement of the nuclei. 

In this paper it is supposed (as was supposed by Heath and Linnett 19480, 
1950) that the orbitals correspond to the hybrids of the central atom, and 
that they are rigid and free to rotate during molecular vibration (Heath 
19486; Cohan 1956). The different modes of vibration are not treated separate- 
ly, however, so that the formulas obtained are invariant under the transforma- 
tions of the molecular symmetry group and are also valid for the combination 
modes. 

THEORY C3, MOLECULES 
The function proposed by Heath and Linnett may be taken as follows: 
V = Vit Voang) 
where 


V, = 4k AP (see Fig. 1), 


A A A 
(1) 2Viane) = Ry [(A’N.A")*+(B’'NB")?4+ (C'NC")*] 


A A A A A A 
+k,[( A’ NB” — A NB)?+(B’ NC” —BNC)?+(C”NA” —CNA)?]. 


Following a series expansion (with w as the dihedral angle between two faces 
of the tetrahedron, see figure) it can be shown that 


A See Je eee Ra hee a cas 
(A'NA")? = : [44"44A""—24 A’. AA" cos (a’—a’")], 


A ae a _ 
(A"NB"—ANB)* = P [4A” cos a’”’+BB" cos (8+w)]’, etc. 
1Manuscript received April 10, 1959. 
Contribution from the Facultad de Ciencias, Exactas y Naturales, Buenos Aires, Argentina. 
Can. J. Phys. Vol. 37 (1959) 
1002 





BLUM AND COHAN: TETRAHEDRAL MOLECULES 1003 





Fic. 1. (a) A, B, C represents the position of H in equilibrium; A’, B’, C’ represents 
the position of H during vibration; NA’, NB’’, NC” represents the position of ‘‘valence line’’ 
direction; N represents the position of N (or similar atom). 

(6) Top view of the figure on the left, inscribed in a sphere of radius / and center N 


are exact to the first degree of approximation. If the variab les are changed: 
A = AA’ cosa’ a = AA" cos a” 
A’ = AA’ sina’ a’ = AA" sina”, etc. 
The formula (1) will be 
(2) 2? Ving = k{A?+4°+B2+B"4+-C4+C"% +0? +0740 +b" +24" 
—2[Aa+A’a’'+Bb+B’b’+Cc+C'c'|+k,[(a+ cos w— 0’ sin w)? 
+(b+c¢ cos w—c’ sin w)?+(c+a cos w—a’ sin w)?} 


and if this formula is differentiated with respect to a,a’, b,b’, c, and c’, the 
necessary condition for a minimum is obtained. Resolving these equations 
leads to: 


(3) 22 Ving = (42+B?+C)Kit+(A"?+B°+C")K2+(AB+BC+CA)K3 
+(A'B’+B’'C'+C'A')Ki+(AA'+BB'+CC')Ks 
+(A’B+B’C+C'A)Ke+(AB’+BC'+CA')Ky 

where 

Ky = [(rit- 2ritritrit-ri—2n41) kypt (hi thithiphol 

Kz = [(l+rit+git2gstritrs—2gs) ket (hi thsthe)ka] 

Ks = 2[(r2—rarstrars—ratst2re—2rire)k y+ (hihot+hohst+hshi)kp] 

Ke = [(rara—rars—2gs+gst2gags—rars)Ry+ (hahs+hshethels)kp] 
Q[(rirstrors—2r3s—regstragstrsgsa—ors) Ret (hihsthshethohs) kp] 


Ks 


Ke = 2[(—rers—rirst2rstrags tres —1sk4—rors) Ret (hog t+hihet+hshs) kp] 
Ky = 2[(rira—2ra—rora—rogstragatragetrors) ke t+ (hohe t+hihsthsh,) Rp] 
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2 


1 , a 2 
es (N+2+x)(N+y); r= D 


(V+y"); 1 = (N+142); 
= 2 (N4+2+x—2'); 


n= Zlte); w= es 5 IVE N (2+) +N (6+2e—2e"—2") 
+y?(2+2x—y"x)]; 


as) i? 
h3 = p & +); 


7; = X (N+2+x—2" ); bk=- 





Nx Ny : 
hy=he= ss = -pW Ey) [N?-+N (8+x—x") + (2+x—2x?—x*)] 
where 
3 > ule 
x = cosw; y = sina; Nz=—; 
kp 


D = N?+N(2+X)+44+2x—2x?. 


To obtain the expression in terms of bending angles of hydrogen atoms, their 
rotation around the central nitrogen atom must be suppressed. Then: 


. Mi... a ee Paz Se 
Fy: = 4 Kit (si+si+53)K2—3¢ Kat (sisets28s+5051)K it @ (si +352—53) 
ial sea : 
XKst+e (—sits2—353)Kote Bsi—srts)K | ; 


: Be seh = a ee 
| = | 2 Kit 2GuetsstssKet dl K+ Gibbs sertasets) 


Sat a ee » 
XKite (sit2ss)Ks+e5 (2sit+s2) Ket, Qsetsi\Ks | 


| l 
$1 = — (x-3); Sp = — (8x—1); $3 = = (l-x 
: by ( Oy ( ) . = ( ). 
TETRAHEDRAL Tz MOLECULES 

The matrix of the potential energy, completely factored, is composed of 
two matrices of 1X1 (corresponding to the modes A; and E) and one of 
2X2 (corresponding to F.). The angular part of the potential function has 
three force constants to be determined, but from the secular determinants it 
is possible to obtain only two constants, even with the data of NDf. This is 
in accordance with the fact that there are only two really independent con- 
stants, which can be initially selected according to the condition of redundance 


A A 
> (R’'NS’)—(RNS) = 
all angles 


Using a treatment completely similar to the one utilized on the NH;-like 
molecules, the following results are obtained: 
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N’ +8 "+21N °"+18N _ 
+10N°+37N 7+60N+3 +36’ 


N + 10N'+34N" +36N 
24(V*+10.N°+37N"+60N +36) ° 


(F.3 denotes the interaction constant between adjacent angles.) 





F,2 = hv 











Fag = k, 


RESULTS 
Infrared frequencies, suitably corrected for anharmonicity, were used in 
the calculation of force constants. The data used are those of Dennison 
(1940) (NH3, ND3), McConaghie and Nielsen (1949, 1953) (PH3, PDs, 
AsH;, AsD;), Haynie and Nielsen (1953) (SbH;, SbD;), and Gamo (1953) 
(NH{, NDf). The geometric constants were taken from the same papers. 
The calculations were performed following Wilson (1955). 


TABLE I 


Force constants 


NH; ND; ° PHs PD; AsHs; AsDs; SbHs " SbD; NH? 








10” X Faz,dynes/em 0.61 0.47 0.72 0.72 0.67 0.67 0.59 0.59 0.2 
10” X Fag —0.129 —0.42 —0.067 —0.068 —0.045 —0.046 —0.035 —0.035 —0.2,—-0.3 
If the values of Table I were put in formula (3), the following expressions 
are obtained: 

For NHs: 
Fog = aqbn(N8+6.351N°+12.586N'+ 15,149? +6.972N° 43.178) 

af NS4-4.858.N°+ 14.280.V'+25.130N°+29.152N?+ 17.635. +5.699 ’ 


k,(1. 131." +3. JOON “+4. 234.V" '+3.834.N" "+1 619N) 


Fas = ¥8 +4.858.V°+ 14.280.V'+25.130N°+29.152.N?+17.635.V +5.699 ° 
For PHs: 

r wba NE+7, LIGN*+ 16.817 N'+23. 294.N*+ 13.889N?+7.354N) 

af N°45.846.V°+20.075.V'+41.358N°+55.607N+37.053.V + 14.644 ’ 

‘ __y(0.01656N°+-0.412.N" ‘+0.710N" '+1.817.N ?40.543.V) 

8 * V°45.846.N°+20.075N'+41.358N°+55.607.V?+37.053.V + 14.644 ° 


For AsH3: 3: 
ky (N°+7.226N +17 A6B7.NV “+25. 105.N° "+15. 157N *+8.023.V ) 


F.: = aah sensqqeet ene eaanicnantnasinaptnauberentaail 


V°+6.006.V°+21.035.V'+44.098.V°+60.162.V?+40.104.V + 16.046 ’ 


h(~ —0.231N °_0.478N ‘0.604. "+0.850N *_0.023.V) 
N®+6.006.V°+21.035.V'+44.098.V°+60.162.V°+40.104.V + 16.046 ° 


For SbHs;: 
k,(N° +7.176N °417.168N “+23.855N" "+14. 569N*+7.7 


Fe = ve tS aes 


~ N®45.932N°+20.589.V'+42.831.N°+58.065.V'+38.729V oat oe 
ky(—0.1198.V °+0.048N *+0.1137N *+1.333N *+0.159.V) 


Fis = 








Fag = N°+5.932N°+20.589.V'!+42.831N°+58.065N?+38.729N +15.424 ° 
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Upon solution of these equations the values of ky and ky, are obtained (see 
Table I1). 











TABLE II 
Orbital following force constants 
NH; PH; AsHs SbH; 
N —8.64 —0.85 —0.32 —0.42 
10” Xk,, dynes/cm 0.89 —1.23 —3.19 —2.01 
10” Xky —7.68 1.041 1.037 0.847 
For the tetrahedral molecules there is only one real solution N = —2 if 


F,2/Fag is negative, and the problem remains, in this case, undetermined 
for ky and ky. 


DISCUSSION OF RESULTS 


It is clear that a negative value of N is anomalous as it means that not a 
minimum, but a maximum of the potential function is obtained; otherwise, 
the type of equation obtained for N allows only negative (real) solutions, 
unless in the case of the ammonia-like molecules, relations between the 
angular force constants are very different from those obtained experimentally. 

The exclusion of a term corresponding to the lone pair of electrons could 
be criticized (Duncan 1953; Coulson 1951), but there are no reliable criteria 
to remove the indetermination of the problem with such a term present. Also 
calculations made by Cohan and Coulson (1956) in which the influence of the 
lone pair was taken into account led to negative force constants. 

On the other hand, in the case of molecules of the ammonium ion type, in 
which these problems do not exist, a single solution for the maximum is 
obtained. 

It is thought that the discrepancy is due to the deformation of the valence 
hybrids during the vibrations. It is possible also, that the valence line approxi- 
mation may have affected the final result. 


LIST OF NON-EVIDENT SYMBOLS 


C3, 174 symmetry group (notation of Mulliken), 
V = potential function, 
V, = part of the p.f. which corresponds to bond stretching, 
Vieng) = angular part of the p.f., 
Al = bond stretching, 
A,, E, Fy correspond to the vibration modes; the notation is the same as the 
; one used by Wilson (1955), 
A” NA’ angle A” NA’, etc., « >" 


A 


a Fd 


AA’ = segment AA’, etc., 
Ll = bond length (see Fig. 1). 
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THE HUND GRAVITATIONAL EQUATIONS AND GEONS'! 


LL. G. CHAMBERS 


ABSTRACT 
Certain solutions of the Hund gravitational equations which represent geons 
are derived, and the stability of steady solutions is discussed. 
1. INTRODUCTION 
In the absence of matter, Hund’s gravitational field equations take the 


form 


(1) vV.g = 0, 
1 og 
tz) vxf = - - at? 
] 9 ‘ 
(3) vt =a (tte), 
- 1 of 
(4) VxXEé = tot fxg. 


By virtue of their non-linearity, these equations admit non-zero solutions 

which are not merely combinations of plane waves, and it has been pointed 

out to me by a referee of a previous paper (Chambers 1959) on this subject 

that such solutions in fact represent geons. It may be shown fairly easily that 

ae conditions on f must apply for a non-trivial solution of equations 
, (2), (3), (4) to be possible. 

— No solution is possible for which f is zero everywhere or for which f 
is the curl of a vector. For were either of these true, the right-hand side of 
equation (3) would be zero, and as it is positive definite, it would follow that 
both f and g would have to vanish. 

(b) Certain conditions on f stronger than those given in the last paragraph 
must be satisfied. For, if equation (3) be integrated over a volume inside a 
bounding surface 


(5) J (f?+g") dr = act f v.tds 


= act f f.dS 
Ss 


and if the last integral vanishes then f and g must be intrinsically zero. This 
will be true for example if S is a large sphere of radius R and the radial com- 
ponent of f is 0(R~*). It may be remarked, however, that equation (5) does 
not preclude the possibility of a geon with a hollow core as for it to apply S 
must enclose v. 


‘Manuscript received April 14, 1959. 
Contribution from the Department of Mathematics, University College of North Wales, 
Bangor, Caernarvonshire, Wales. 
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2. STATIC GEONS 


Static solutions of equations (1) to (4) are given by fs, $s, which satisfy 


(6) V .gs = 0, 

hs 
(8) V.fs = 52 (fstgs), 
(9) VX£8s = + fsXs. 


In general, it is not possible to solve these equations. If, however, either (a) 
f; and gs are parallel or (8) $s vanishes, a solution is possible. In both these 
cases equation (9) becomes 


(10) Vv Xg;5 = 0. 
Thus from equation (7) 

(11) fs = 2°Vu 

and from equations (6) and (10) 

(12) @s = 2?Vv 

where 

(13) V4 = 0, u and v are dimensionless. 


Now, in case (a), fs and gs are parallel and this implies that Vu and Vv 
are parallel. That is « and v form the same equipotential surfaces. In other 


words 


(14) u = F(v), 
(15) fs = 2c°F'(v) Vv, 
(16) V .fs = 2c?F''(v)(Vv)?*, 


using equation (13). Substituting in equation (8), 
2c? F"’(v) (Vv)? = 2c*{ (F’(v))? +1} (Vv)? 


whence 
F’(v) = tan (v — v9) 


where vp is a constant, 


whence 
(17a) f; = 2c? tan(v — v)Vv 
(17d) gs = 2c°Vv 


where v is subject to equation (13). On differentiation vo drops out and 
v—vU is also a solution of equation (13), and so the solutions may be written 








i 
p 
ff 
r 
f 
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(18a) fs; = 2c? tanvVv = 2c*¥ (log sec v) 
(185) és = 2Vv. 


This method indicates how the case (8) may be dealt with. Putting 
fs = 2c?F'(v)Vv = 2V{ F(v)}, 
where v again satisfies equation (13), equation (8) becomes 
2c? F''(v) (Vv)? = 2c7{ (F’(v))*} (Vv)? 


whence 


F'(v) = 





Yo—v- 
Here again the constant v may be dropped, and 
2c" 2. 
(19) f,= “a Vv = —2cV(log v). 
As the range of solutions of equation (13) is large a correspondingly large 
set of solutions is possible for fs and gs, subject always to the conditions of 


Section 1 being satisfied. 
The equivalent mass distributions are given by 


1 2, 2 1 
20 ~-——, s) = -7—V-fs. 
(20) ree (fst+gs oy * Ss 
In case (a) the equivalent distribution is 
(21) ciate V. {2c* tanv¥v} = ae sec’ v(¥v)” 
4ry 2ry 


. a . . . . . 
and in case (8) the equivalent distribution is 


1 = \ c (vv)? 
99 So 6 aig fa ee 
fee) 4rry = v = 2ry vu 
Using Green’s theorem the total equivalent mass is 
1 f 1 f 
2 can eee ft 4 1 SS . 
( 3) dry ly s) dr dry « a dS 


where S is the bounding surface. 


3. PERTURBATION OF HUND’S EQUATIONS 
Hund’s equations are 
(24) vV.g =0, 


25) vVxf = —-~, 
. 


Bi 
(26) V.f = —4ryptsza (fits), 
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ss 1 0f 4zxypv, 1 

9 ss ys a ; 
(27) vx c al +3 (fxg) 
For a steady universe, equations (25) and (27) become 
(28) ¥ Xf _ 0, 

: 4rypv , 1 
(29) VX8s = ——_— +23 (fsX8s). 


p and v are independent of time. 
If f and g be disturbed, p and v remaining the same, the disturbance may 
be discussed in terms of F and G where 


f = f,+F, ¢ = £s+G. 
Substituting these in equations (24) to (27) it follows that 


(30) Vv .gs+V.G =0, 


1 aG 
(31) VXfst+VXF = oa 
(82) Vist ¥.F = —dryptya { (fst F)'+(8stG)'}, 
8) VX BSF UXG = GAT A Ist F)X Bot G)I, 


Subtracting out the stationary parts from equations (30) to (33) it follows 
that 


(34) vV.G=0, 
5 a 
(35) VXF = at 
(36) V-F = +) [fs FH8s-G] +5 (F' +6"), 


1aF 


I 
~ gp tz (F X8stfsxG+FXG). 


(37) vVxG = 
It will be observed that equations (34) to (37) do not depend explicitly on p 
and v, but only implicitly through fs and gs. Thus the discussion which 
follows will be true both for geons and for mass distributions. 

In order to discuss the stability of f and g, it is necessary to linearize 
equations (36) and (37). These then become 


(38) v.F = “: (f5.F +85-G) 
and 
(39) wxG = 22E 14 (Pxest fsxG). 
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Equations (34) and (35) are solved immediately by putting 


(40) G = VXA, 
0A 

(41) F = ms 

where 

(42) ct = 7. 


(It is not necessary to introduce a scalar potential (Stratton 1941a).) 
Writing 
asc’ = fs, Bsc? = $s 


equations (38) and (39) become 


a 0A 
(43) <a (vV.A) = as. 5 +6s- (VXA), 
(44) VX(VXA) = oh sx eA asx (WXA). 


Equation (44) may be simplified slightly. From equation (28) 


(45) fs = CVu 

and it follows that if 

(46) i a asi 

equation (44) may be rewritten as 

(47) vxtW(w XA} = —Woh+1g.x24. 


Extending the ideas of non-linear mechanics (Minorsky 1947) it is possible 
to consider each stationary configuration of f and g as a point in an abstract 
phase space. When there is a disturbance, a path is traced out in phase space 
according as to the type of stability of the system at that point, and the 
temporal variation in the linearized problem. 

The possible types of steady fields will be considered and the conditions 
found which A must satisfy. 


4. NATURE OF STABILITY OF STEADY SOLUTIONS 


(a) Stable Vortex Point 
In this A = Ay exp(twr) (Apo real). Equation (48) becomes 


(48) iw{V .Ao—as. Ao} = 6s.(V X Ao), 
(49) Vv KF W(V¥ XAv)} w*WAot+iwW Bs X Av. 


I 


Splitting into real and imaginary parts 


(90a) Vv a Ay— Qs. Ay — 0, 
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(506) Bs.(V XAo) = 0, 
(50c) BsXAo = 0, 
(50d) Vv X{W(¥ XAv)} = w?WAo. 


It is also necessary that w* is positive. Multiplying equation (50d) scalarly 
by Ao and integrating over the whole of space it follows that 


(51) w fa (dr = fac. [Vx {W(¥xX Ao)}] dr 


J Wwx Ao)? dr— f W[AoX ( XAv)]-dS 


and if this last integral vanishes, w? is always positive. If this is not the case, 
then the splitting up into real and imaginary parts is not admissible. It will 
be noted that w vanishes with VX Apo and if this be the case all the equations 
(50) are satisfied identically. This case is not stable, however, as it corresponds 
to a time variation Ct+D. 

It may be remarked that (50c) sometimes implies (500) for 


0 = V. {8s X Ao} = Ao.(VX Bs) —Bs.(VX Ao) = —85.(V X Ao) 
by virtue of equation (9) if fs and gs are parallel. 


(b) Focus 
In this A = Age??e™7 (Ao real) where p is positive for instability and 
negative for stability. 
Equation (47) becomes 
(52) VX {W(VXAo)} = (b+tw)*WAot+(p+iw) Ws X Ao. 
Taking the imaginary part of equation (52), it follows that 


(53) 2pwAot+iwBsX Ao == (). 


This is a vector equation the only possible solution of which is Ay = 0. Thus 
there is no steady state which corresponds to a focus. 


(c) Saddle Point 

In this case the time variation is either e*”*? depending on the nature of 
the small disturbance. 
(d) Stable Node 

In this case the time variation is e~:” or e~?:? where p; and ps2 are positive, 
depending on the nature of the small disturbance. 


(e) Unstable Node 

In this case the time variation is e?:7 or e?:7 where p; and pz are positive, 
depending on the nature of the small disturbance. 

These three cases, (c), (d), (e), may be dealt with by assuming a time 
variation e*? in equations (43) and (47). 








1014 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


Equation (43) becomes 
(54) k(V.Ao—as.Ao)+6s.(VXAo) = 0 
and equation (47) becomes 
(55) VX{W(¥ XAo)} = —k?>WAl+RWBsX Ao. 
Multiplying scalarly by Ao 
(56) Ao. (VX{W(V XAv)}] = —RWAS. 
It follows therefore that there cannot be nodes, because the two roots of (56) 
are equal and opposite in sign, and therefore there must be a saddle point. 
Thus the stability of a steady-state solution depends on the quantity 
(57) Ao.[VX{W(VXAo)}] = W(VXAo)?—V. { WAoX(VXAp)}. 


If this is positive there is stability and a vortex point. If it is negative, there 
is a saddle point and instability. 
It follows therefore that a condition for stability is 


(58) {W(¥ XAo)2dr > [W[AoX (VX Av)].dS 


and in particular this is true when the integral on the right-hand side vanishes. 
This may be done by any of W, Ao, VYXAo vanishing on the surface, by the 
surface being parallel to either Ay or VXAo, by the product WAyX(¥ X Ao) 
being sufficiently small at infinity if the whole of space is under consideration, 
or of course by some orthogonality property. 

Further conditions for stability which must hold are equations (50a) to 
(50c) in the general case, equations (50a) and (50c) when fs is parallel to g5 
and equation (50a) when $s vanishes. 


5. AN EXAMPLE OF A GEON AND ITS STABILITY 
For simplicity consider a static geon in which $s vanishes, and there is 
spherical symmetry. Consider the potential 


a 
59 one 
or ) ¢ . 
( r 
This function is harmonic and 
° 7 2 
te, a0 2e°. 
(60) ia xg Oe 
v r 


The value of W corresponding to this is given by 


(61) WV = e-*, f; = Vu. 
Clearly u = —logv? 

and so 

(62) W = r*/a?. 
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The equations governing a perturbation A are given by 


9 
(63) v.A—" 4, = 0, 


9 . 
(64) VX (VXA)—= rX(V XA)—w A = 0. 
These equations can be solved by the following method (Stratton 19410). 
Let 
(65) A = YX(Lr) 
where L is a scalar. It is clear that this equation satisfies equation (65) 
identically as 


- o es: este en _ _1ob 
(66) A,=%, Mere, 6h 


Expanding equation (64) into its components it is found that the 7 component 
vanishes, and that the 6 and ¢ components are respectively given by 


diac ia dl i L (: cine Eye om |S =) 
mn rsin 6 t Oolr sin 6 on sin @ dg or-od 


Fists S). 6 
v 


r ror sin 6 d¢@ rsin@d¢ 
and 
(68) thee. a 1 (a + oy 
r \ar'a0 a6 Lr’sin 8 (2 Gne ae ag/ If 
214 .) w” aL 
“yr a (- 0) + 0° 
It will be seen that equations (67) and (68) are satisfied provided that 
aL 2aL i ae aL 1 aL 
" mL 1 @ = 
(69) art r arte sin 6 00 (ein @- 00 a? sin 0 Oe te! a 0, 
i.e. VL+u"L = 0; 


a solution of this equation is given by 


(70) L = {AghS” (wr) +A oh (wr)} Sinn (8, &) 


where h{ and h{” are spherical Bessel functions (Stratton 1941c) and S,,,(@, @) 


is a surface harmonic. 
The stability depends upon the value of the surface integral in equation (51). 


If the surface is a sphere of radius R, 


4 
(71) fwiax(wxa)}.d8 = *, fiaxcwxay], ae 
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where dQ is an elementary solid angle. This becomes, on substituting the 
appropriate values for A and VXA 


" ef | a (ay ( 1 ab)" | 
2 peek a ace gs oar a eee . 
Ae) a” aa dr \\ a6 7 sin 8 0¢ = 
It is clear that in general this does not vanish. 
However, if L = Lojn(wr) P,(cos 8) where Lo is a suitable dimensional factor, 


expression (72) becomes 


Tne: ’ ain Sed 
(73) 2 * Lowjn, (wR) jn (wR) fair, (cos 8) sin 6} 
and if either j,(wR) or j,(wR) vanishes expression (73) vanishes and the 


disturbance will be stable. 
Thus in general the geon considered will be unstable, but will, however, 


be stable for certain perturbations. 
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ON ‘‘PLASTIC’”’ WAVES IN SOLIDS! 


JoHN M. BowsHER? 


ABSTRACT 


The study of the propagation of ‘‘plastic’’ waves in solids has reached a stage 
where it is necessary to consider which direction future research should take. In 
the past 90 or so years many experiments, mostly designed to elucidate certain 
points of engineering significance, and a few attempts at a theoretical study have 
cast some light on the subject and revealed it as one of formidable difficulty. 

Nearly all the experiments have of necessity relied on rather dubious theories 
for their interpretation, and part of the present paper will be devoted to a descrip- 
tion of an apparatus which gives results capable of being interpreted with a very 
minimum of theory. The remainder of the paper is devoted to a short review of 
past work with particular emphasis on basic phenomena and to a brief discussion 
on the most pressing problems still remaining. The experiments described in the 
present paper bring to light a factor in the propagation of ‘‘plastic’’ waves that 
seems to have been overlooked in previous work. 


INTRODUCTION 

Most work on the subject of the propagation of stress waves in solids has 
been devoted to the limiting cases of what shall be called here Low Amplitude 
and Very High Amplitude waves. In low amplitude waves, the stress amplitude 
is assumed to be much less than the elastic limit so that, after the passage of 
the wave, the specimen practically regains its original form. In other words, 
low amplitude waves are merely elastic or sound waves. When the stress 
amplitude is so great that the conventional elastic constants cease altogether 
to have any meaning and the specimen material flows like a fluid, the wave 
motion is referred to as very high amplitude. 

In High Amplitude stress wave propagation, however, the stress amplitude 
is comparable with (say 1 to 100 times) the elastic limit, and is not so large 
that the detailed mechanisms governing the flow of solids can be ignored as 
with very high amplitude waves. The chief characteristic of propagation in the 
high amplitude range is the very large amplitude dispersion the input signal 
suffers. Frequency dispersion has long been familiar in the low amplitude range 
where the concepts of group and phase velocity have been thoroughly explored 
and meaningful values of the physical constants of the specimen may now be 
given after measurements of the group velocities only. In the high amplitude 
range nothing has been found in the literature on frequency dispersion and this 
is, therefore, a problem which should be tackled soon. Amplitude dispersion is 
not so well known since it is an effect of very small importance in low amplitude 
work, but in high amplitude waves the velocity of a stress increment situated 
above the yield stress of the specimen material may, in certain cases, be as 
little as 10°% of the elastic wave velocity (Wood 1952). 

An input signal consisting of a step-function-like stress pulse thus rapidly 
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loses its initial shape and becomes transformed according to the properties of 
the propagating medium and the amplitude of the initial pulse. For conven- 
ience, we shall use the word profile to describe the shape assumed by an 
initially step-function-like pulse. Mathematically, the profile may be con- 
sidered as the o’, u, €, or p versus C curve, where 


o’ is the engineering stress 
u the particle velocity 

e the strain 

p the density 

C the stress wave velocity, 


and the present paper is devoted entirely to consideration of the profiles of 
high amplitude pulses. 

Previous work on the profiles of high amplitude waves is wholly theoretical 
and is typified by the figures (1, 2, 3) to be found in von Karman and Duwez 





Plastic wave front 





Elastic wave front 
4 


fc, a 


Fic. 1. Typical theoretical profile of a stress pulse propagating in a wire (after von Karman 
and Duwez 1950, used with permission). 


(1950), Taylor (1946), and Broberg (1956). The profiles by von Karman and 
Taylor are for a medium having a stress-strain curve always concave down- 
wards, and, moreover, having no strain—-time effects. Broberg’s profile applies 
to a medium exhibiting the delayed yield phenomenon and also applies to 
stresses large enough for the stress-strain curve to be concave upwards leading 
to the development of shock fronts. 

Since so many assumptions had to be made to derive these profiles it 
appeared valuable to design an apparatus capable of obtaining directly the 
profile of a stress pulse propagating through a medium. The apparatus, which 
is to be described later, obtains the u—C profile directly but it has proved 
possible (Appendix) to derive some simple, but approximate, correction terms 
to the well-known relation from elasticity, ¢ = pCu, where a is the stress. The 
corrected equation may then be applied to a high amplitude pulse propagating 
in a medium with an unknown stress-strain curve. The corrections are based 
on the Lagrangian formulation of ‘‘plastic’’ waves by von Karman but are 
independent of strain—time effects since all measurements are reduced to those 
capable of being made while the pulse is actually being propagated. 

The particle velocity u is measured by assuming that the free surface 
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Fic. 3. Typical theoretical profile of a stress pulse propagating in a medium exhibiting a 
stress-dependent time delay for yielding (after Broberg 1956, used with permission). 


velocity after a normal reflection of the pulse is 2u, which assumption has 
been shown to be valid within very small error up to stresses much in excess 
of those used in the present apparatus (Walsh and Christian 1955). 
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A NOTE ON PLASTIC AND “PLASTIC” WAVES 
Since nothing is yet known about the front separating the elastic and plastic 
regions in the medium through which the wave is propagating, the term 
‘‘plastic’’ wave is used here throughout. It is felt that the term plastic wave is 
inappropriate as the region is not plastic but is an intermediate stage between 
elastic and plastic. 


LOW AMPLITUDE WAVES 

It was estimated recently by Davies that well over 1000 papers have been 
published since World War II on the subject of stress waves in solids; obviously 
it has become well nigh impossible to survey the field in any detail, and this 
section is intended to be just an introduction to give the reader some references 
where he may find further bibliography. We shall discuss briefly papers on the 
particular problem of the impact of two bars only since results obtained from 
experiments on impacting bars give information on some of the phenomena of 
importance in high amplitude stress wave propagation. 

The first papers, by Cauchy (1826) and Poisson (1833), are not of much 
value at present and it was undoubtedly de St-Venant (1867) who gave the 
first satisfactory solution of the problem of the impact of bars with plane ends. 
Simultaneous impact over a plane end is virtually impossible to realize in 
practice, however (Schneebeli 1871; Hamburger 1885), and Sears (1908, 1912) 
in his work used a combination of de St-Venant’s theory and Hertz’s (1882) 
formula for the deformation resulting from the compressive contact of two 
surfaces with radii of curvature R; and Re. Hertz himself did briefly consider 
the problem of impact in the appendix to his paper and obtained results similar 
to those of Davies (1949) but that work seems to have been forgotten until 
recently. 

The experimental arrangement of two bars, one being hit by the other, is 
extremely suited to the study of short, large, stress pulses and considerable 
understanding of plastic flow has been gained by using the Hopkinson Bar, as 
a refinement of the two-bar apparatus is known. Hopkinson (1914) estimated 
the time integral of a pressure pulse by measuring the momentum transferred 
to a detachable “‘time-piece”’ at the end of a pressure bar, but recently Davies 
(1948) has applied modern techniques to the bar and an illustration of their 
use may be found in Taylor and Davies (1942), where the dynamic yield stress 
of Cordite was measured. Flow of the pressure bar material takes place quite 
frequently in Hopkinson bar experiments. 















HIGH AMPLITUDE STRESS WAVES 


The subject of high amplitude wave propagation is intimately bound up 
with that of plasticity and yield phenomena, and a favorite method of 
investigating plastic phenomena is to impact the specimen and study the 
resulting “‘plastic’’ waves and/or flow. Any study of ‘‘plastic’’ waves must, 
therefore, be preceded by some study of plastic flow and the various effects 
found. These will not be discussed in detail but the references quoted give 
typical results and also lead the reader to a greater bibliography. 
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Let us first list which effects are of importance in the study of “‘plastic’’ 
waves: (7) strain rate, (iz) temperature, (777) pressure, (iv) delayed yield, 
(v) attenuation. 


(t) Strain Rate 

For many years it has been realized that for many metals, particularly iron 
and steel, a different set of mechanical constants applies when the strain rate 
is high, and very many experimental investigations have been performed to 
elucidate the difference. As early as 1872, Hopkinson (1872a and 6), the father 
of the inventor of the Hopkinson Bar, found that an iron wire could withstand 
a momentary stress of almost double the elastic limit without any noticeable 
impairment of elastic properties and since then there have been many experi- 
ments devoted to attempts to evaluate the differences between static and 
dynamic properties. The best general reference is Rinehart and Pearson (1954) 
and a typical result for iron has been found by Whiffin (1948) to be 
D/S = 3 for S = 20 tons/in*?, and D/S = 1 for S = 120 tons/in* where D is 
the dynamic yield stress and S the static; these results show that the problem 
is of considerable complexity. Many attempts have been made to explain the 
strain-rate effect on the basis of an “‘incubation”’ period for yield but they 
overlook a paper by Edwards, Phillips, and Jones (1940) in which it is reported 
that for iron at very low strain rates, the yield stress begins to rise again. 

The strain-rate phenomenon is not only important in metals; in plastics also 
it is found that the response to a stress depends on the rate of loading, which is 
important now with the steadily increasing use of rubbers and high polymers 
for vibration isolation. The physical properties of high polymers with long 
chain molecules with, perhaps, various degrees of cross-linking may, fortu- 
nately, be represented by a relatively simple ‘‘memory” function and may now 
be calculated with some confidence (cf. Alfrey 1948 for example). 


(it) Temperature 

Increasing strain rate and decreasing temperature have qualitatively 
similar effects on many mechanical properties of metals but no simple relation 
has yet been found between increase of strain rate and decrease of temperature 
since every metal and every alloy has a different behavior. Attempts have been 
made by Holloman (1946), however, to obtain a single-valued function that 
will express the behavior of a metal under various conditions of strain, strain 
rate, and temperature but his equations have only a limited application. 

A fairly comprehensive survey of the effects of low temperatures on the 
mechanical properties of metals has been given by Russell (1931), who gives 
the following summary. 

With decrease of temperature for common metals, the 

vield point increases, 
yield strength increases, 
reduction of area decreases, 
impact resistance decreases, 
hardness increases, 
and the endurance limit increases; 


the close similarity with rapid loading effects may clearly be seen. 
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There are, though, many exceptions to the above summary, notably the 
stainless steel alloys (Krivobok 1952), which exhibit greater ductility at low 
temperatures. 

The effect of temperature on the mechanical properties of a metal must 
therefore be established for each metal and also as far as possible under the 
conditions for which knowledge of the property is desired, since present knowl- 
edge is inadequate to allow a reliable prediction from data obtained under one 
set of conditions to another. 


(i111) Pressure Effects 

The effect of pressure on the elastic constants of solids is not very great in 
the range of pressures of present concern. Hughes and Kelly (1953) and Cook 
(1957) give examples of the use of Murnaghan’s (1937) theory to the calcula- 
tion of the velocity of low amplitude waves in a solid under hydrostatic stress, 
and the reader is referred to those papers for further details. 


(iv) Delayed Yield 

With some materials, there is an appreciable time delay before yield takes 
place when the specimen is loaded to a stress above its yield point. This delay 
depends on the conditions and may be as much as 6 or more seconds and is of 
great importance in impact problems since the mechanical response of the 
material to a transient load is going to depend very strongly on the time con- 
stants associated with that load. 

A great deal of both experimental and theoretical work has been carried out 
on the delayed yield phenomenon, chiefly because it occurs to such a large 
extent in those materials, such as iron and steel, which are of great use in 
engineering. The basic theoretical paper is that by Cottrell and Bilby (1949) 
which, together with annotations by Yokobori (1952a and b), describes a 
mechanism of yielding in terms of dislocations in the crystal lattice of iron and 
dislocation atmospheres round each iron atom formed by the disturbing effect 
of the smaller carbon atoms which tend to cluster in the low-potential regions 
around the iron atoms. In this theory a time delay in yielding is predicted. 
Typical experimental results are quoted by Wood and Clark (1949, 1951) and 
Vreeland, Wood, and Clark (1953), who measured the time of yielding by the 
number of repeated short-duration pulses required before yielding took place. 
The 1951 Wood and Clark investigation found that the effect of temperature 
upon time delay for mild steel confirmed the theory of Cottrell and Bilby. 

The properties of crystals which are likely to be of importance in high 
amplitude pulse propagation and which have not been discussed above are 
fully dealt with by Cottrell (1953), whose work includes a comprehensive 
bibliography. 


(v) Attenuation 

For some reason the fact that the ‘‘plastic’’ wave will suffer attenuation 
during propagation and that this will affect conditions at the front between 
the plastic and elastic regions quite markedly seems to have escaped attention 
until now, and it is believed, therefore, that the discussion following is the 
first in print. 
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PROPAGATION OF HIGH AMPLITUDE WAVES 

The present theories of the propagation of high amplitude waves such as 
those by Taylor, von Karman, and Rakhmatulin are treated in standard works 
(for example, Kolsky 1953), and will not be discussed here save to emphasize 
that they are only extensions of ordinary elasticity theory to finite strains; 
there has not yet appeared in the literature a rigorous definition of a plastic 
wave and a subsequent analysis of its behavior for a given set of conditions. A 
physical interpretation of plastic waves presents grave conceptual difficulties 
which, we feel, are very inadequately considered in the literature and for this 
reason all references here to high amplitude waves are to “‘plastic’’ waves and 
not the loose term plastic waves which can lead to confusing and possibly 
inaccurate assumptions on the modes of displacement taking place during 
the passage of the wave. 

The three theories mentioned above apply only in cases where the stress— 
strain curve is always concave downwards. When the curve may become 
concave upwards further difficulties arise due to the formation of shock fronts 
(see, for example, White and Griffis 1947, 1948). 

If the strain-rate phenomenon is included the problem becomes still more 
complicated and quite complex mathematics have to be employed to obtain 
solutions to the simplest case (Alter and Curtis 1956; Malvern 1951; Rubin 
1954). A profile calculated by Alter and Curtis is shown in Fig. 4; note that ina 


Mox. Slope 
we yy? 
of NS 






Strain 


Wx/a 1 x/a, to*x/o, +7 
-_-_oOo-->> 
Time-t 


Fic. 4. Typical theoretical profile of a stress pulse propagating in a medium exhibiting the 
strain-rate effect. The dotted curve shows the result of a two-step application of the same 
pulse amplitude (after Alter and Curtis 1956, used with permission). 


material exhibiting the strain-rate effect plastic strain components can travel 
with the velocity of sound, Co, in contrast to the rate-independent model. 

Experimental investigations into ‘‘plastic’” waves are very few and are 
rather too strongly dependent on a particular theory to be entirely satisfying 
(see Kolsky 1953). A puzzling paper by McReynolds (1949) describes a type 
of “‘plastic’”’ wave unlike the usual type considered here; it is difficult to explain 
McReynold’s waves but they may be an experimental artifact. 

In the very high amplitude range the stress amplitude is so great that the 
specimen material may be treated as a fluid with a considerable simplification 
of the theory. It is not uncommon for measurements to be made on two 
parameters of the motion simultaneously in very high amplitude work (Walsh 
and Christian 1955; Doran, Fowles, and Peterson 1958), which enables the 
equations of the motion to be solved without having to assume an equation 
of state. 
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In high amplitude work, on the other hand (Duwez 1947; Duwez and 
Clark 1947), no experiments on propagation involving the simultaneous 
measurement of two parameters had been reported when this work was begun 
(but see Bell 1959). It seemed, therefore, that an attempt should be made to 
record the profile of a high amplitude stress pulse in a way which would be 
independent of a particular theory, and could be said to mean something by 
itself. To do this, two parameters at least must be measured simultaneously. 


APPARATUS 


Let us first list the techniques that are available to one wishing to measure 
high amplitude waves. The problem divides naturally into two parts: (a) 
generating devices, and (0b) receiving devices. 


(a) Generating Devices 

The only really convenient method to generate a high amplitude stress wave 
in a solid is to hit the solid with a suitably designed hammer. The use of high 
explosives is unnecessary because the stress amplitudes will then usually be far 
too large. Methods based on using the pressure developed by the passage of 
an electric spark through sea water contained in a small volume are rather too 
inflexible for experimental work where a wide range of pressures is desired. 

The design of a hammer depends on many circumstances, some of which 
have no scientific significance. In the present work, for example, the commonest 
method of driving a hammer (allowing it to fall freely under gravity) was not 
available for the simple reason that the laboratory ceiling was not high enough. 
The device actually used was a high-speed pneumatic hammer (Fearon 1953), 
as this was thought to be a more flexible and reproducible tool than a catapult 
in which the hammer is accelerated by rubber bands. The hammer used was 
shaped with a spherical end of 5-ft radius of curvature, large enough to ensure 
essentially plane waves over the measuring area. A flat-ended hammer cannot 
be used because the alignment difficulties become too acute (vide supra). 


(b) Receiving Devices 

Most techniques are ill-suited to the measurement of two variables at once, 
but some worth considering are: 

1. Walsh and Christian's method (1955).—This will not work in its original 
form in the high amplitude field, but the basis of the method is capable of 
adaptation (see later). 

2. Pins.—This is the conventional technique for very high amplitude work 
but is not suitable for high amplitude work because the setting-up tolerances 
become too fine. Another disadvantage is that the apparatus is destroyed after 
each experiment. 

3. Piezoelectric transducers: (A) accelerometer type, (B) pressure-sensitive 
type.—The main difficulty with these devices lies in their calibration. They 
are known to become increasingly non-linear as the stress amplitude increases, 
but as yet no one has found a way to calibrate over the whole range of ampli- 
tudes since at the top of the range the crystal fractures and is lost. Minshall 
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(1955) has reported the use of Z-cut tourmaline crystals for very high amplitude 
work and it appears that they would be usable in the high amplitude range. 
Only one parameter is measured, though, a(t). 

4. Resistance strain gauges.—These are too large for the sort of measurements 
at present of interest. 

5. Diffraction gratings ——This extremely elegant method has been described 
by Bell (1956). The main disadvantage is that a lathe accurate to optical 
standards is required to prepare the specimens. 

6. Condenser microphones.—The condenser microphone is immediately 
attractive as there are no problems of bonding the specimen to the microphone 
(as in 3); it is merely necessary to have a conducting plate a short distance 
away from and parallel to the specimen surface. The use of a condenser micro- 
phone to record the surface motion of a specimen is well known but another 
parameter was needed. Obviously, one cannot measure the stress or density 
directly with a condenser microphone, so a way must be found of measuring 
the velocity of the stress pulse itself as it moves through the specimen. 

Following the method of Walsh and Christian a technique involving the use 
of a small inclined groove cut in the surface of the specimen was devised for 
the present work; this is the first time any such method has been attempted. 

Thus the receiving device used consisted of two condenser microphones 
one measuring the surface displacement of the specimen, and the other the 
stress wave velocity in the body of the specimen by means of an inclined 





groove. 

Calibration was performed in the following way. 

(a) The surface motion (W)* system. The associated electronic circuitry 
was designed so that the response extended to direct current and the micro- 
phone was calibrated directly by comparing the amplitude of the movements 
of the surface after the passage of a stress wave with movements of the micro- 
phone plate measured by a micrometer. Thus the calibration is independent of 
linearity or distortions in the recording apparatus (both electronic and optical) 
and is referred directly to a micrometer screw. 

(6) The stress wave velocity (C) system is calibrated in an entirely different 
way. It is necessary to get some clear indication of the time the stress wave 
reaches various points in its traverse of the specimen. This was done in the 
following manner: let there be on the conducting surface of the microphone 
some regions extending transversely right across the surface where the con- 
ducting layer is removed and let the microphone surface be inclined so that it 
is parallel to the face of the groove. As the stress pulse reaches the groove, the 
voltage on the microphone will start to change due to the change in capaci- 
tance produced by the reduced plate-specimen spacing and will continue to 
change until the ‘‘front’’ of the pulse has reached the first break in the con- 
ducting surface. When the front of the pulse reaches the next effective region, 
there will be a change in the working area of the microphone and the output 
voltage will have a sudden change in slope. Thus if the breaks in the conduct- 
ing surface are at accurately known intervals and the times between the 


*The letter used for displacement in the present work is w. 
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discontinuities on the record can be measured, the velocity of the stress pulse 
can be determined directly. 

A disadvantage of this arrangement is that for accuracy a very large band- 
width is required. The microphone used was made with two breaks in the 
surface about 3 mm apart (see Fig. 5), and the groove and microphone were 





Fic. 5. Diagram showing the arrangement of the two condenser microphones under the 
surface of a specimen. 


inclined to the surface of the specimen at an angle of about 10°. Thus the two 
breaks are roughly 0.4 mm apart in the direction of travel and since a typical 
longitudinal stress pulse velocity is approximately 5 mm/ysec the separation 
in time of the breaks in the recorded trace is about 0.1 usec. To obtain 1% 
accuracy of measurement, therefore, a system rise time of about 1 nanosecond 
is needed. It was not possible to construct such a system but the method was 
thought to be sufficiently promising to be worth developing and a trial 
apparatus was constructed with a time resolution of about 30 nanoseconds. 
Thus the C system measurements made are subject to considerable error but 
are mainly of interest in establishing the feasibility of the method. 

The rise time of the W system was under | usec and was low enough to be 


of no concern. 


EFFECTS ON “PLASTIC”? WAVE PROPAGATION OF THE VARIOUS 
PHENOMENA DISCUSSED ABOVE 

Various profiles may be predicted qualitatively by assuming reasonable 
values for the parameters and then drawing out successive phases of the motion. 
It is virtually impossible to predict the profiles quantitatively because so little 
is known about such factors as attenuation, strain rate, velocity of unloading 
waves, frequency dispersion of high amplitude waves, etc. Another difficulty 
is that some materials exhibit a polymorphic transition at certain stress 
amplitudes; these transitions are well known in both static and very high 
amplitude work and there is no reason to assume their absence in high 
amplitude work ab initio. A series of experiments on ice may, perhaps, be of 
value in a study of the effects of polymorphic transitions on the propagation 
of high amplitude waves. 


| 
| 
| 
| 
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The author (Bowsher 1958) has drawn some profiles qualitatively and 
initial results on lead show that some of his predictions have been confirmed. 
In particular, it was found that when the effects of attenuation are considered 
for an amplitude dispersive medium (as all media are above the yield point) 
the elastic front preceding the ‘‘plastic’’ one is ‘‘maintained”’ by the latter at 
the yield point and suffers no apparent attenuation until the “plastic” front 
has decayed to a value only slightly above the yield point. This comes about 
in the following way: for a ‘‘plastic’’ front to propagate in a stable manner, the 
specimen material immediately in front of it must be prestressed to the yield 
point. In a medium where the stress-strain curve is composed of two straight 
lines joining at an angle at the yield stress (case 01000 of Bowsher) the 
“plastic” front travels at a fixed velocity (<Cpo) and, therefore, can only 
maintain the necessary prestress by ‘‘sending” out small elastic waves every 
time the attenuation has reduced the amplitude of the elastic front sufficiently. 
An unloading wave is sent back towards the impact point every time this 
happens to maintain momentum balance. Thus the profile of the elastic 
precursor is changed because the regions behind the front are maintained at 
stresses above those which would pertain if simple attenuation were the only 
mechanism effective. Similarly, the region behind the ‘‘plastic’’ front is at a 
lower stress than simple attenuation would predict. It appears on the profile 
that the “plastic’’ front is maintaining the elastic one. 

When attenuation and dispersion are present (i.e. the o’—e curve is curved 
above the yield point, case 01100 of Bowsher), “‘plastic’’ waves of small 
amplitude can travel at velocities approaching Cy and the “‘plastic’”’ front can 
maintain a stable condition without having to generate elastic-maintaining 
waves to precede it. The velocity of the ‘‘plastic’”’ front gradually increases as 
the wave becomes attenuated. 

When an attenuating medium also exhibits the strain-rate effect as formu- 
lated by Malvern, plastic strains can travel at Co and the elastic front is 
maintained at the yield point (case 11000 of Bowsher). In his experiments the 
author was not able to obtain a “plastic’’ front with a velocity greater than 
about 2Co/3, which indicates that strain rate has much less effect than 
dispersion in his experimental arrangement and suggests that the reciprocal of 
the strain-rate parameter for the material investigated is somewhat less than 
the value suggested by Malvern of 1 usec, and may be as small as 0.1 usec. The 
attenuation measurements show, though, that the strain-rate effect is not 
entirely absent since the elastic front is maintained at the yield stress 
(see Fig. 9). 


EXPERIMENTAL RESULTS 


Some results have already been mentioned above in the discussion on 


attenuation phenomena of ‘‘plastic’’ waves. 

The value of an apparatus that can measure two parameters simultaneously 
is shown by the profiles in Figs. 6 and 7; the latter, the u—C profile, is entirely 
experimentally obtained and is not dependent on the assumptions of any 
theory (other than that the free surface velocity is 21). 
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Fic. 6. Experimental o—C profile for lead. 
Fic. 7. Same results as for Fig. 6, plotted as a u-C protile. 
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Fic. 8. Typical experimental w-t curve. The point marked ‘plastic wave arrival’ shows the 
arrival at the free surface of the elastic wave reflected from the “‘plastic”’ front. 

Fic. 9. Attenuation curves for elastic and ‘‘plastic’’ waves in lead. The upper lines are for 
‘‘plastic’’ waves, and the lower for elastic. Points marked @ are for measurements made at 
+20° C and points marked O for measurements at —30° C. 
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The C system measurements are not of very great accuracy since, as men- 
tioned above, the equipment was designed merely to test the feasibility of such 
an experimental approach. The results show that the method can be used 
provided a short rise time for the impact is achieved—this is chiefly a matter 
of experimental technique. The accuracy of the system depends almost 
entirely on the resolving time of the electronic circuitry following the con- 
denser microphone and may be very great indeed if the latest developments 
in subnanosecond rise time amplifiers are used. The accuracy obtained in the 


present experiment is only +15% but this is sufficient to show the worth of 


the method. 

The W system was found to be capable of measuring u to within less than 
3% error which is sufficient for the determinations of profiles in the present 
work. On the experimental w-t curve shown (Fig. 8) the errors are too small for 
inclusion on the graph, but are: w +1 yw, ¢ +0.1 usec. 

The profile shown (Fig. 6) was obtained for lead specimens cast from a sheet 
of ‘“‘chemical’’ lead and it may be seen that the yield stress is 97-410 bars, and 
the elastic wave velocity 1.25+0.1 mm/usec. 

It is not meaningful to give attenuation results numerically as the effects 
discussed above tend to obscure the simple attenuation phenomenon. The 
results obtained (Fig. 9) show clearly that the elastic front is indeed maintained 
at the yield stress until the ‘‘plastic’’ front is no longer able to do so. This 
occurs at a specimen thickness of 20 mm for the lead and impact velocity used. 
It is interesting that the only difference found between experiments performed 
at +20° C and those at —30° C is expressed in the attenuation results. In all 
other respects the specimens behaved in a closely similar way at the two 
temperatures; it will be valuable to extend the scope of the temperature effect 


measurements to elucidate this point. 


CONCLUSIONS 

A discussion has been presented on the propagation of “‘plastic’’ waves in 
solid specimens. The main intent has been to summarize the main effects for 
the benefit of other workers in the field and to attempt a study of the physical 
aspects of “‘plastic’’ wave propagation. 

An apparatus has been described which obtains results independent of 
theoretical assumptions and helps one establish which of the many phenomena 
are likely to be of importance for any particular specimen material. 

The work presented here shows that there are many factors of ‘‘plastic”’ 
wave propagation about which virtually nothing is known experimentally. 
Among these are: 

1. Velocity of unloading waves in a solid. This is always assumed to be the 
same as that of elastic waves (Co). There is no a priori reason for this assump- 
tion though static stress-strain measurements do indicate that the velocity 
may be C». 

2. The detailed microphenomena in a 


‘plastic”’ front. 


3. Attenuation. The author’s experiments are exploratory in nature but do 
seem to confirm his mechanism in which the amplitude of the elastic precursor 
is maintained at the expense of the “plastic’’ front amplitude. 
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4. Temperature effects. The author’s experiments are once again only a 
beginning and a detailed study would be most valuable. 
5. A confirmatory series of experiments on other profiles, e.g. p—C or e—C, 


would be of interest. 

6. Frequency dispersion of high amplitude waves. Nothing is reported in the 
literature on this and in the author’s experiments the specimen length was 
deliberately kept short (<30 mm) to minimize any frequency dispersion and 
so no information was gained on this topic. 

7. The present experiments have only been on lead; other solid materials 
should, of course, be investigated. It would be of particular interest if experi- 
ments on a material like solid argon were performed since this is one of the few 
solids whose lattice structure is very simple and amenable to theoretical 


calculation. 
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APPENDIX 
DERIVATION OF A RELATION BETWEEN STRESS AND PARTICLE VELOCITY IN AN 
ELASTIC-PLASTIC MEDIUM 
Let us consider a material having a stress-strain relation like that shown in 


Fig. A. Let the elastic limit be at P with co-ordinates (0, ¢.) and let the input 
stress have amplitude o, resulting in a strain €,. Note that the o’~e curve is the 





e Pp 


Fic. A. Stress-strain curve for an arbitrary medium. 


one applying in the dynamic state of the material under the influence of the 
pulse, and also that, in general, it is different to the one used by von Karman 
since we shall be assuming an infinite extent of medium. 

In Lagrangian co-ordinates the particle velocity u is given by 


‘Vd do’ 


ee i Y¥oa* 
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where p is the density of the medium (assumed constant). This relation, strictly, 
only applies at the end of the medium where it is impacted but to the present 
accuracy it is sufficient to assume that it applies throughout the extent of 
travel of the pulse. 
Let 7(e) be the slope of the region between P and Q, 
M _be the slope of the region between O and P; 


then 
ue= | Vil/o de 
0 
oe na ap ea 
and Up =| Vil/o det | VT/p de. 


€ 


Now o = f “pC? act foci de, 
where pCi = M; pCi =T; Ci = Cie), 
and Uy = fa at fc: de 
s e 
= Cut J Ci de, 


a‘p 
o> = pCicet J, pCi de. 
a‘p : 
Thus pCwy = pC Cree t+ pCi J. C, de 


’ Cu a’p a*p 

Tp ) . 2 

and 2 * av = cuter ~eorog | Ci de-Cif Ci ae| ° 
p *e £e 

If the region of the stress-strain curve from P to Q is linear, the term in square 

brackets becomes 


Ci(ep— €e) _ Ci(ep— €e }= 0, 


i.e. the first term in the right-hand side is the correction term to allow for the 
change in slope at P. Therefore the second term is the additional correction 
factor allowing for the curvature from P to Q. 

To get an approximate evaluation of the term in square brackets let us 
assume DeJuhasz’s stress-strain law. We have already assumed p = constant 
so that e must remain fairly small. 

We assume 
for Oe Sa, oa’ Me 


ll 


ll 


Ge S &, oa’ = BvVe 
eet Crees o, = Me = BV e&, 


therefore, e. = B*/M?. 
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Now, the velocity of ‘‘plastic’”’ waves is 





1 do’ 
Os p de 
= a 4/2. aa 

2pv/e 2p‘ 

Thus [ | becomes 
“p —7i ep ate 

J, ae a-4/2 in /2cau- 4/& [ engl et de 
substituting for B, 
y 2 2 2 
i.e. { j= Cied ig-a-2 ats a]. 
Now € = Ue/Co 


€ 4 ‘ i 
therefore { j= ciut| « i(1 -; 4)+(4) \2(#) -1}| ; 


€, is given by a) = BY 6; 


92 72 92 


p Dp 
€p = See re meee 
B €ep Co Uep Co 


' : 2 a te\'S2(ue\? ,\ 
therefore [ ] = Cou} Es ~utican)+(%) B a =i |. 


Thus 





eee at, 2 ae a(#) 
ed ong i(- 2 tira) - pC ity + pute (Co— C)+oCan {2(#) if 


Uep 


3 = 2/3 1/6 2/3 
= ee Chas *! (p—Ue) 


SUMMARY 
Thus for an elastic wave we have 
Te = pColle, 
for a “‘plastic’’wave with a stress-strain curve consisting of two linear portions 
(1) = = pCotte+pCi(up— Ue), 


whilst for a parabolic stress-strain curve 


3\** aya. 16 2/: 
= Cw+(3 pCi'*Chue (tp — te) e 


bo 
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To compare these formulas we substitute some figures given by DeJuhasz— 
these may not necessarily be correct but they give some idea of the magnitudes 
of the corrections to be expected. 


Let us take p = 8.9 g/cm’, 
Co = 4.5 mm/usec, 
C, = 0.9 mm/usec, 
ue = 100 cm/sec, 
uy = 150 cm/sec. 


Then the terms become 
pCoute = 401 bars, 
pCi(up—uUe) = 40 bars, 


3 2/3 . . 
(3 pCi! *Chur!® (4p —ue)*”* = 20.5 bars. 


Thus eq. (1) gives 441 bars, and eq. (2) gives 421.5 bars, a difference of 5%. 

A recent paper by Feltham and Meakin (1957) leads one to suppose that the 
linear relation, eq. (1), might be more appropriate for copper in the present 
stress range, and we use it throughout. 

When the stress amplitude is large enough for a shock wave to develop the 
Rankine—Hugoniot equations apply, and it is found that stress and particle 
velocity for a “‘plastic’’ wave are related by the formula 


(3) oy = pCity. 
Let us rewrite eq. (1) as 
o) = pCityp+ptte(Co— Ci). 

It may thus be seen that eq. (3) is a limiting case of eq. (1) for u, > u, and 
C,; — Cy the conditions applying in a shock wave. For example, using typical 
figures we have u, = 10%, C,—Co = Co/10. 
pCiu, 

100 ’ 


i.e. only a 1% change is introduced by our correction terms. 


Whence op = pCiuy+ 


NOTE ADDED IN PROOF.—An interesting mathematical treatment of this problem by L. W. 
Morland (1959. Phil. Trans. Roy. Soc. London, A, 251 (997), 341) has just appeared. 
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ABSTRACT 

Rotational energy levels of axially asymmetric nuclei have been calculated in 
the manner of Davydov and Filippov, and a comparison with experiment shows 
good quantitative agreement. It is suggested that agreement may be improved 
by introducing the Bohr—Mottelson vibration-rotation interaction and a centri- 
fugal stretching correction analogous to the type used in molecular spectra. The 
D-F method seems to be particularly useful for nuclei in the transition regions 
between “rotational’’ and ‘‘near-harmonic’’ modes of collective excitation. 


Bohr and Mottelson (1953) have shown that even-even nuclei with sphe- 
roidal shapes may be expected to exhibit rotational states with energy levels 


given by 
Exot = FL I([+1 I=0,; j 
“rot ~ 23 ( + De ( ’ a 4, 6, 


where & is the effective moment of inertia and J the total nuclear angular 
momentum. The ratios of the energies represented in this simple level scheme 
are, E(4+)/E(2+) = 3.33, E(6+)/E(2+) = 7.00, E(8+)/E(2+) = 12.00. 

Davydov and Filippov (1958) have extended the treatment to include 
nuclei with ellipsoidal shapes, i.e. nuclei not possessing axial symmetry. 
Deviations from axial symmetry may be characterized by the nuclear shape 
parameter y, a quantity which varies between 0 and 2/3. Values of y of 0 
and m3 characterize prolate and oblate spheroids respectively. The Hamil- 
tonian operator for rotational motion is given by 

i? 3 i 
he ao 2 2 sin’ |y— (27/3)i} 

where B is the mass transport parameter, 6 is the nuclear quadrupole deforma- 
tion parameter, and the I; are operators of the projections of the total angular 
momentum along the body-fixed axes of the nucleus. The Hamiltonian is 
symmetric about y = 7/6 and consequently a prolate ellipsoid of deformation 
y (<a/6) will exhibit the same rotational spectrum as an oblate ellipsoid of 
deformation 7/3—y 

The energy equations of Dennison (1931) for an asymmetric top may be 
applied. The symmetry conditions of Bohr (1952) exclude certain energy 
levels. Resultant curves of relative energy vs. y are given in Figs. | and 2 
The 2+, 4+, 6+, 8+,... spin sequence appears and, as long as y is small, 
the relative spacing of these levels does not differ appreciably from that of a 
spheroidal rotor. New levels, of energy E(2’+), E(3+), E(4’+),..., appear, 


'Manuscript received June 29, 1959. 
Contribution from the Radiation Laboratory, McGill University, Montreal, Que. 
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Fic. 1. Rotational energy levels of axially asymmetric even-even nuclei. 
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Y DEGREES 
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Fic. 3. Variation of the nuclear shape parameter y~ with mass number 4. Figures indicate 
proton number. 


however, as a direct consequence of axial asymmetry. If, as proposed by 
Davydov, y is determined from the ratio E(2’+)/E(2+) and the experi- 
mental value of E(2+) assumed, the energies of the remaining levels may be 
predicted for even—even nuclei with large quadrupole deformations. A com- 
parison of these predictions with experimental data in the range 150 < 4 
< 250 is shown in Table I. All spin assignments are those given by the authors 
of the reference quoted. In general, agreement is good—especially in the 
tungsten—osmium transition region and in the region of heavy nuclei. In 
these regions discrepancies are typically of the order 1 to 3%, with the 
exception of the 4’+ level of Os!® where the discrepancy is 13%. 

Figure 3 shows the variation of y with mass number A for the range 
150 < A < 250. Since the 2’+ level is not always known, y has been calcu- 
lated from the ratio E(4+)/E(2+) in some cases. y calculated in this manner 
is always somewhat larger than that given by the E(2’+)/E(2+) ratio. The 
three transition regions from ‘‘rotational’”’ to ‘‘near-harmonic”’ level schemes 
(Scharff-Goldhaber et al. 1958) are clearly shown for A ~ 150, 190, and 
220. The samarium-gadolinium group (A ~ 150) exhibits a very fast transi- 
tion, whereas the tungsten—osmium group (A ~ 190) and the radon-radium- 
thorium group show a more gradual transition. 

Qualitative arguments indicate that the minimum finite value of y for a 
stable deformation cannot be small, since the frequencies of vibration and 
rotation would then be comparable, and separation into rotational and 
vibrational levels meaningless. The data in Fig. 3 would seem to support 
this. However, most of the low values of y were obtained from the E(4+)/ 
E(2+) ratio. It is well to keep in mind that low values of y obtained from 
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ratios of the energies E(2+), E(4+), E(6+), E(8+) do not necessarily prove 
the existence of a stable deformation of this magnitude. For values of y up 
to approximately 15°, the effect of assuming a y is equivalent to introducing 
the Bohr—Mottelson vibration-rotation interaction correction of the form 
—bI?(I+1)?. Values of y obtained from the E(2’+)/E(2+) ratio are not 
subject to this objection. 

Table II gives the results of calculations of energies of what has been 
termed the ‘‘ground state rotational band’’, 2+, 4+, 6+,8+,.... If y is 
determined not from the ratio E(2’+)/E(2+) but rather from any two 
members of the “ground state band’’, then the energies of the remaining 
two levels of this band may be predicted to an accuracy usually well within 
1%. However, energies of the 2’+, 3+, 4’+,... levels cannot then be pre- 
dicted with any precision. The data is consistent with assuming a higher 
value of y for the ‘‘ground state band” than for the 2’+, 3+, 4’+,... levels. 

A qualitative argument may help to clarify this point. Two corrections 
which may be of importance in the Davydov—Filippov treatment are the 
vibration-rotation interaction correction and a centrifugal stretching cor- 
rection. As in molecular spectra these corrections have the form —b/?([+1)? 
and become more important as the equilibrium deformation decreases. For 
the ‘ground state band” it is the 6 deformation that is important and for the 
2’/+,3+, 4’/+,... levels it is the y deformation that is important. The 
effect of this correction to the energies of the levels of the “ground state 
band”’ is equivalent to an upward shift of y of several degrees. Since y 
is relatively insensitive to the ratio E(2’+)/E(2+), the correction which 
must be applied to the 2’+ level does not appreciably affect y as determined 
from this level. Thus the difference between y calculated from the “ground 
state band” energies and y calculated from the E(2’+)/E(2+) ratio is quite 
naturally accounted for by the type of correction arrived at by analogy with 
molecular spectra. 

This paper has shown that the Davydov and Filippov treatment can 
successfully account for a large number of nuclear energy levels. It cannot, 
however, predict all the low-lying excited levels, for it is known that other 
0+ and 2+ levels exist which cannot arise as excited states of a rigid asym- 
metric rotor. This point has been clearly shown in a recent paper by King 
and Johns (1959) concerning the energy levels of Os!**. 
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IRREVERSIBLE THERMODYNAMICS AND KINETIC THEORY 
IN THE DERIVATION OF THERMOELECTRIC RELATIONS! 


J. W. LEECH? 


ABSTRACT 


The variety of experimental behavior found in thermoelectricity at low 
temperatures provides an incentive for extending the theoretical treatment. 
In the past, irreversible thermodynamics has been used to provide a satisfactory 
derivation of the Thomson relations. It appears that what one might calla general 
‘primary equation” of thermoelectricity can be derived from irreversible 
thermodynamics, and that this should be of direct use in kinetic discussions. 
The form of the equation suggests that an even more direct approach to the 
thermoelectric problem should be possible. 


1. INTRODUCTION 


The phenomenon of thermoelectricity has for a long time furnished a 
convenient method for temperature measurement. It is less well known as a 
means for investigating the solid state though recent experimental investi- 
gations (see for example MacDonald et al. 1958a, b and 1959; cf. also Mac- 
Donald 1959) have revealed a variety of new phenomena, particularly at low 
temperatures. There have been two main theoretical approaches to the 
subject. On the one hand, irreversible thermodynamics has been used (see 
for example Onsager 1931; Callen 1948; de Groot 1951; Domenicali 1954) 
to provide a satisfactory derivation of the fundamental Thomson relations; 
on the other hand, kinetic theory based on the Boltzmann equation has been 
used for detailed calculation of the various thermoelectric quantities (see 
for example Wilson 1953). If, however, we combine these approaches there 
results an appreciable simplification of the kinetic calculations. The generality 
of the irreversible thermodynamic argument suggests that it may always be 
profitable to use it as a check in the development of the kinetic explanation 
of the more complicated thermoelectric phenomena. 

We believe it worth while to summarize both treatments using a common 
notation before combining them to obtain a simplified expression for the 
calculation of the Thomson coefficient. This should obviate the rather surpris- 
ing misidentifications and unphysical statements which are to be found in the 
literature of the subject* and which can usually be resolved by an under- 
standing of both forms of argument. 

1Manuscript received May 13, 1959. 


Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada. 

Issued as N.R.C. No. 5300. 

2Now at the Department of Physics, Queen Mary College, University of London, London, 
England. 

*Mott and Jones (1936, p. 309), for instance, appear to identify the Peltier heat with the 
contact potential difference between two media and ter Haar (1956) writes down expressions 
for the electric and heat currents (his equations (38) and (39)) which depend on the thermo- 
electric power in an obviously wrong manner. 


Can. J. Phys. Vol. 37 (1959) 


L044 







LEECH: IRREVERSIBLE THERMODYNAMICS 


2. IRREVERSIBLE THERMODYNAMICAL TREATMENT 


We wish to take into account the simultaneous flow of heat and electricity 
through a medium which, in general, contains both thermal and electric 
potential gradients. It is supposed that, although the medium is in a non- 
equilibrium state, it is nevertheless possible to define regions small on a 
macroscopic scale but large on a microscopic one whose instantaneous states 
may be specified using thermodynamic co-ordinates. 

Irreversible thermodynamics can be approached from a number of different 
directions and at a number of levels (see de Groot 1951; Denbigh 1951; 
Prigogine 1955; Davies 1957). Entropy plays a dominant part in all these 
approaches, particularly its creation and transport. We propose to intro- 
duce the irreversible thermodynamic postulates in the form of an equation 
linking the entropy flow J‘) and the entropy density s by means of the 
balance equation: 

















5 (s) I 
(1) a= Ved teh a. Me 





Here, the J, are flow vectors (in our case the electric and heat currents) 
and the X, are thermodynamic ‘‘forces” giving rise to these flows. The term 






1 
> JX; 





obviously represents the rate of creation of entropy per unit volume. 
It is further assumed that, for conditions not too far removed from equili- 







brium, we may write: 











(2) J, = ~ 1,sX; 


J 







where the /,;; are functions of temperature and of the constitution of the 
medium which obey the Onsager relations: 










(3) lig = 14. 









Assumptions (1) to (3) appear in most treatments of the subject though in 
many cases they are deduced from more fundamental postulates. In parti- 
cular the Onsager relations (3) are usually held to be based ultimately on the 
principle of detailed balancing. 

There is considerable choice in the precise form of either the flow vectors 
J; or of the appropriate conjugate quantities X;. However, once one set has 
been chosen there is no longer any arbitrariness in the other. 

In the case of thermoelectricity the choice of the electric current density 
J and the heat current density Q as flow vectors is clearly indicated since 
they play a major role in the kinetic treatment. The corresponding conjugate 
quantities, the thermodynamic forces X,, will be identified by setting up an 
entropy balance equation based on first- and second-law considerations. 
Restricting attention to one-dimensional flow, the time variation of the 
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internal energy density of the electrons in a conducting medium is, from the 
first law of thermodynamics: 


(4) oe = ~ 82 (1,8) 


where ¢ is the electrostatic potential at the point in question. 
The second law leads us to write (see Sommerfeld 1956): 











(5) 2 i 


where g is the chemical potential per unit conduction charge in the neutral 
medium, e the electronic charge (in sign and magnitude), and N is the 
electron density. 

It will be seen that consideration has been restricted to a single type of 
charge carrier but that no specific mechanism has been assumed for the 
process of heat conduction. 

Combining (4) and (5) and allowing for the continuity of charge we find 
after some rearrangement: 


© fo —mltolra@)le-rgen 


where ©, = — 0¢/ 0x is the electrostatic field acting on the electrons. 
The following identifications may now be made: 


(7) the entropy flow JS = - {0.—gJe}, 


the thermodynamic ‘‘force’’ conjugate to J; isX y = ©,—T 2 (g/T), 
(8) : 


see ” . Re a Out a 
the thermodynamic force conjugate to Q,isXg = T— BG) ‘ 


hence, assuming near equilibrium: 


ee et ae’ 
Jr = lyy & T= (g/1 ) +100] 7 3 i. 
6. = tes 72/7) [400 2. (4)] 
r Qu z Ox g, QQ ax \T 


where the Onsager relations (3) show that: 


(9) 


(10) lag = Iz. 


This completes the statement of the postulates of irreversible thermo- 
dynamics in a form relevant to our considerations. It remains to identify the 
various thermoelectric quantities and to determine their necessary inter- 
relations. 


LEECH: IRREVERSIBLE THERMODYNAMICS 
3. IDENTIFICATION OF THERMOELECTRIC QUANTITIES AND 
RELATIONSHIPS 


Equations (9) may be rewritten in a form more appropriate to the later 
kinetic discussion as: 


J; 


a’g 
aT 
‘a. dx 
x 
where the symbol 0’g/ dx denotes change of g at constant temperature, i.e. 
it allows for the possibility of a non-homogeneous medium. In terms of the 
newly defined coefficients the aaa reciprocal relation (10) becomes: 


\ 


(12) p= oT 4T ZR e/D]. 


We also note that the thermal conductivity is given by: 


sia : | a7 CE-y) 


and that for a homogeneous medium: 


(14) = [J:/Ez]arjaz—0, 


which justifies our use of this particular symbol. 

The identification of the Thomson and Peltier coefficients is made by con- 
sidering dq/ dt, the rate of heat absorption by the electrons per unit volume. 
Under steady-state conditions J, = const. and dq/ dt = du/ dt; hence from 


(4): 
(15) 


For a homogeneous medium (0’g/ dx = 0) we get using (11) and (15): 


: ag _ Se he | cle rt. 
(16) a J: 7 97 (8/2) |+ = 


The Thomson coefficient » is defined as the reversible heat absorbed per 
unit volume when unit electric current flows down unit temperature gradient. 
This immediately gives: 


(17) - a4 (B/c) 
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from the realization that reversing the electric current or the temperature 
distribution leaves the first and third terms of (16) unchanged but alters the 
sign of the second. This is essentially the identification made by Mott and 
Jones (1936) though they made it by assuming the term d/dx [x(d7T/dx)] to 
be negligible. 

The identification is in fact a precise one and involves no assumptions beyond 
those basic to irreversible thermodynamics. 

Equation (17) has been used in the kinetic approach to calculate the Thom- 
son coefficient. In that form no explicit use has been made of the Onsager 
relation. If this is now remedied by combining (12) with (17) we get: 


(18) w= To) 7S 


This result does not seem to have appeared previously. It represents a useful 
practical advance on the form (17) since the quantity 8 has been replaced 
essentially by g. Whereas in the kinetic approach the determination of 8 
requires the evaluation of a separate integral, g is to be identified with the 
Fermi energy and can be assumed known. 

The Peltier coefficient 4, is defined as the reversible heat absorbed at the 
junction between two media A and B when unit electric current passes from 
A to B. Our formulation cannot be used to deal with a strictly discontinuous 
boundary but we may consider the more realistic case of a continuous merging 
of medium A into B under isothermal conditions. For such a region, again 
assuming a steady state, we obtain using (11) and (15): 

dq J; d 


c foe oy es 
(19) ars Jz 7. (8/0 g). 


Or, integrating from an interior point of A to an interior point of B: 


cy fern fmallt-LE-a]} 


The division into reversible and irreversible parts again occurs (this time 
solely in terms of current reversibility) and the Peltier coefficient is identifiable 


as: 
(21) Han = [8/o—gla—[8/o—g]p 
which, using the Onsager relation (12), becomes: 
is a E ae. | E a6. 
9 SERS So) Sees 
(22) r oar B ot ar as 
Combining (18) and (22) now gives the familiar Thomson relation: 
d | Bas Ma—~ HB 
9) amas) saree ae See 
(23) arvl or rT 


From the derivation it is clear that the Peltier heat is associated with the 
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transfer of charge from the interior of one medium to the interior of the other 
and is not, as sometimes asserted, directly connected with surface charges or 






contact potentials. 

Finally it is necessary to obtain expressions for the thermoelectric power. 
To do this we must first evaluate the open-circuit potential difference between 
two points such as X and Y of Fig. 1, which represents the usual arrangement 








T >T, 






Fic. 1. Basic thermoelectric circuit. The region including the contacts X, Y below the 
dashed line is supposed to be isothermal, as are also the corresponding contact regions 
at temperatures 7) (supposed fixed) and 7. 







for measuring thermal electromotive forces. This is done by integrating the 
electric field between X and Y along a path through the media A and B. The 
temperature 7; is supposed fixed in what follows while 7 may be varied. 
The integration is effected by again considering boundary regions as con- 
tinuous. The transition region from the interior to the exterior of medium A 
at X produces an effect equal and opposite to that at Y since they are both 
uniformly at the temperature. This leaves integrations over the temperature 
range 7; to 7 in media A and B and over the transition regions from A to 
B at temperatures 7; and T. We define the thermoelectric force of the couple 











as 








Eg—E, = Ey—Ex in Fig. 1 (and we note that 7 > 7}) 





= -f €,dx. 





Using (11) with the restriction J, = 0 we find after some working: 





oT 
(24) Eg—E, = Han(7)—Han(Ti) + f (un—wa) AT. 
T, 





The thermoelectric power of the couple is given by 












(25) Se— Sx = (Ex— Ex) 






and we can then show using (23) and (24) that 


_(dSs dS 
(26) t (43s ay = (up—Ha) 





and that 







IIp—U, = 7(Sg—Sgp), writing Ip, = Mg—TIq. 





It is then obviously consistent with these relations to define an absolute 
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thermoelectric power S, and Peltier heat, II, for a single medium related to 
the Thomson coefficient, u, by the equations 


dS 
7 ae = 
(27) = Var’ 
and 
ican, 


This completes the derivation of the Thomson relations. It should be 
emphasized again that all the results are independent of special assumptions 
regarding the boundary regions of the media and rest upon distinguishing 
between reversible and irreversible heat production. The somewhat artificial 
nature of the latter distinction should also be noted. For observational reasons 
J, and dT/dx are regarded as independent variables in the above argument 
but any two of J,, dT/dx, €,, and 0; could, in principle, be so chosen and 
each pair would lead to a different division of the heat production into rever- 
sible and irreversible parts. 


4. THE KINETIC APPROACH 

The simplest calculations of the magnitudes of thermoelectric quantities 
are based on a kinetic approach which starts from the Boltzmann equation 
and takes into account only electron transport processes (see Wilson 1953, 
Chapter 8). 

We suppose it possible to divide the medium into elementary regions which 
contain sufficient electrons to be characterized by a distribution function 
I(x, ¥, 2, Pr, Py, Pz) giving the distribution of electrons in phase space. This 
function then obeys the Boltzmann equation: 


od af eres (2) 
(28) ant .vf+p. V,f = at) con’ 


The following additional assumptions are now made: 
(a) there is local equilibrium, i.e. 


(b) there exists a unique relaxation time 7 such that: 


at) _ _f=fo 
sa (2 coll ~ - 


where fo is the distribution function corresponding to over-all equilibrium 
when no gradients or flows are present; 

(c) the only forces acting on the electrons other than random collision 
forces are those due to the electric field, i.e. 


(31) Pp; = eG. 
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With these and assuming a variation in f in the x direction only, the Boltzmann 
equation may be rewritten: 









of 


(32) f = fo—r, L_ Pee 


m Ov, ° 






The electric and heat current densities may be written in terms of f as: 






Je = e fofa’s 





(33) 





. 
OQ, J Linv’ v, f d’v. 


Using (32) and (33) in conjunction with the properties: 







(34) 








electron density, 


we find: 






(35) 


’. . . m 20 
0, = -F G, Pee v ae -> fre of d’y. 






It is difficult to determine f explicitly, but approximate values for the 
integrals in (35) can be obtained by assuming f = fo, the equilibrium value 
of the distribution function. A second approximation is necessary in this 
determination. Since fy is known only as a function of T it is required to assume 
a value for the temperature at the point in question before evaluating the 
integrals. However, since true thermodynamic equilibrium does not exist there 
is no strict meaning to temperature. The difficulty may be overcome by 
assuming a temperature equal to that which would exist in the same body 
in true thermodynamic equilibrium if it had the same mean value for some 
such property as the kinetic energy of the electrons. With this sort of assump- 












tion (35) becomes: 






aT 
Ao er op Sth. dx 





(36) 





, tf, 2 oat g mt [irytp gt OF 
a9 Urt an, 2 9 & 2 oT vz fodv. 7. 









In this form contact is made with the irreversible formulation since equa- 
tions (36) are of the form (11) with (0’g/ dx) = 0. The irreversible method is, 
of course, more general than the kinetic one though the main conclusions still 
hold and, in particular, it is possible to identify the Thomson coefficient as in 
(18): 
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oT CS Or: 
with 
Ne’, 
¢= cae 
(37) 
“ i ee 
a= er fof fod’ = an NCa 


where Cy is the electronic specific heat per electron, i.e. 


a= Te S26, der 
aT 6 aT) * 


The thermoelectric power for the medium is similarly: 
2C., @ 
(39) - (26u1_ 28), 


5. PARTICULAR CASES 


Two elementary cases may be considered: 


(38) 


(a) Metallic conductors.—For these g is the Fermi energy of the electrons 
per unit positive charge; hence, g = ¢/e where ¢ is the Fermi energy per 
electron. 


Now, 
2, 27.2 
ark 
§ £0 12 f 
where 
h? ay" 
= 8m \ x 
its value at 7 = 0, and 
w’k?7 
Ce “2r0 : 
therefore, 
es wkT : _ Cer 
(40) S= De ta and aie 


(b) Semiconductors.—Here 
2 = ‘ a\—3/2 
¢=kT In 15 h’(2% mkT) 


and C,,; = 3k/2, assuming a Maxwellian electron distribution, hence: 


Sk ¢ *. . 
(41) S= (t2) and = = “ 5 


a 
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MacDonald (1958) has previously shown that the result un = C.y/e follows 
quite generally from the kinetic argument with the condition of eq. 30 but 
without specific assumptions concerning the electron statistics. The order of 
magnitude of the result is correct over a wide range of conditions both for 
metals and semiconductors, showing that the assumptions introduced in the 
kinetic treatment are reasonably valid. The relevance of Thomson's remark 
(Thomson 1854) to the effect that ~ may be called the specific heat of 
electricity is readily apparent. 

So far the electron transport process has been supposed unaffected by the 
crystal lattice through which the electrons move. It has been suggested 
(Gurevich 1945, 1946) that the electrons may receive directed momentum 
from the lattice waves. The effect on the thermoelectric power may be seen 
by a small modification of the above simple kinetic argument. 

Let us, following MacDonald (1954), assume that the lattice waves exert a 
pressure, p, on the conduction electrons proportional to uy, the energy density 
associated with the lattice waves, i.e. p = euz, where e is of order unity. The 
net force density is therefore given by 















a dT 
ar (eur) 7 





for a homogeneous medium. Assuming a fraction \ of this force acts on the 
electrons then the accelerating force per electron becomes: 


X98 (ay,) 2, 
N ar "Gx ° 


substituting this in the Boltzmann equation leaves the coefficient ¢ unchanged 
but modifies a to: 


(43) + a a Gis N aT é (a,)]. 


On the assumption that no net increase in the internal energy of the electrons 
results from the electron/lattice wave interaction the irreversible thermody- 
namic argument remains unchanged and the thermoelectric power becomes: 


_ S21 * >}. # 
(44) S=\3'e a = 













(42) F, = eG€,— 









mat 





6. CONCLUSIONS 










It has been shown how irreversible thermodynamic considerations simplify 
the determination of the various thermoelectric coefficients using kinetic 
theory arguments based on the Boltzmann equation. In particular an ex- 
pression for the Thomson coefficient is derived (equation (18)) which does 
not seem to have appeared previously in explicit form. Although only the 
simplest models have been considered, the extension to more complex ones is 
reasonably straightforward. It is immediately apparent, for instance, from 
equation (11) that the same formal kinetic relationships must hold for a 
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non-homogeneous material as for a homogeneous one if ©, is replaced by 
{€,—(0’g/dx)}. The restriction to a single type of charge carrier can also 
readily be removed. 
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NOTES 


THE ENERGIES AND RELATIVE PAIR PRODUCTION CROSS SECTIONS 
FOR Zn® AND Na?? GAMMA RAYS 


P. P. Stincu, H. W. Dosso, AND G. M. GRIFFITHS 


For gamma rays just above 1.02 Mev a three-crystal pair spectrometer of 
even moderate resolution can be used to determine the gamma-ray energies 
very accurately. Since exactly the rest mass of the pair electrons is subtracted 
from each incident photon by pair production, a measurement of the small 
remaining kinetic energy of the pair electrons establishes the energy of the 
gamma rays with considerably greater accuracy than that of the kinetic 
energy measurement. This method has been used to determine the energies of 
the gamma rays from Zn® and Na* sources. The spectrometer consists of three 
sodium iodide (ThI) crystals as shown in Fig. 1. The annihilation of the 


GAMMA RAYS 


Single Channe! Single Channe! 
Analyzer Anolyzer 


KICK SORTER 


Fic. 1. Block diagram of the three-crystal spectrometer. 


positron after the pair event in the center crystal (1.75 cm X 4cm X 4 cm) 
gives two .51-Mev quanta which leave the center crystal in opposite directions. 
The colinear .51-Mev annihilation quanta are detected in coincidence in the 
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side crystals (2 in. long and 1} in. diameter). The pulses produced in the side 
channels enter single-channel differential discriminators which select only those 
pulses which correspond to .51-Mev absorption in the side crystals. The out- 
put from the side-channel analyzers and the pulses from the center channel 
pass into a triple coincidence circuit. The triple coincidence output pulses 
open a gate which allows coincident center channel pulses to enter the 30- 
channel kicksorter. The differential discrimination in the side channels greatly 
reduces the background due to coincidences from multiple scattering, double 
Compton events, and cascade gamma rays and so increases the resolution of 
the apparatus. 

The three Nal (ThI) crystals were mounted on RCA 6342 photomultipliers 
and had a resolution of about 8.38% for 1.3-Mev gamma rays. A 0.2-in. 
diameter collimated beam of radiation, obtained with a 3l-cm long lead 
collimator, was incident along the center crystai axis. The collimation 
improved the resolution significantly. Photomultiplier gain changes were kept 
to a minimum by using stabilized power supplies and relatively low counting 
rates. No significant gainshifts were recorded over a period of 1 day. The 
kicksorter spectrum was recorded after every 2 to 3 hours and thus a regular 
check was kept on the gainshifts. Linearity of the center crystal was established 
to better than 0.1% for the energy range from 40 kev to 500 kev by the use of 
a standard pulse generator and gamma-ray sources of known energies from 
Cs87, Eu, and Na”. 

Typical spectra produced by Co® alone, and by Co® with Na”, and by Co® 
with Zn® are shown in Fig. 2. Since the Co® y-rays have been established to 
better than 1 kev (1.8325+.0003 and 1.1728+.0005 Mev), using a double- 
focusing spectrometer (Lindstrom 1953) to compare the Co® lines with the 
accurately known RaC 1.4158+.0002 Mev transition, the energies of the Na” 
and Zn® y-rays can be obtained from the present experiment to better 
than 2 kev. 

The individual peaks for each of the gamma rays were isolated graphically 
by subtracting the typical Co® spectrum and the background from the com- 
posite spectra as shown in Fig. 2. Assuming a Gaussian shape for the peaks, 
the mean energy and the variance of each peak was computed with the help 
of the University of British Columbia Alwac IIIE Computer. The statistical 
uncertainty in the mean energy for each peak was also calculated from the 
variance and total number of counts in the peak; in all cases this uncertainty 
was not greater than | kev. A more realistic estimate of the uncertainty in the 
mean gamma-ray energies was obtained by calculating the standard deviation 
in the measured means from seven runs for the Na® and Zn® sources. Since the 
Na” and Zn® gamma-ray energies were determined by simultaneous measure- 
ment with respect to the Co® gamma rays, gainshifts and other systematic 
errors were effectively cancelled out. The final results of the present experiment 
are 


Zn® gamma-ray energy 1.1124+.0019 Mev, 
Na” gamma-ray energy 1.2736+.0018 Mev. 









NOTES 





o-0-o- (0) 


eeveee (b) 


(c) 








COUNTS 
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KILOVOLTS 





Fic. 2. Typical spectra obtained with three-crystal pair spectrometer. Curve (c) is for 
Co, curve (6) for Co® and Zn®, and curve (a) for Co® and Na?* gamma rays. 






The present results for Zn® gamma rays agree, to within the errors, with 
those reported by Johanson (1956), 1.112+.003 Mev, and Waggoner (1950), 
1.112+.007 Mev; and they are significantly different from those of Hedgran 
(1950), 1.125 Mev, and Good (1951), 1.127+.009. The Na® gamma-ray 
energy agrees with that of P. Marmier as reported by Ajzenberg (1955), 
1.275+.005. In both cases our estimated errors are less than those reported 









earlier. 









Relative Pair Production Cross Sections 

We have also measured the relative pair production cross sections for y-rays 
of Zn®, Na*, Co®, and RdTh, using the pair spectrometer. For a source of 
cascade gamma rays such as Co* the relative pair production cross section is 
simply obtained from the ratio of the number of counts in the two pair peaks 
since the number of gamma rays of each energy is the same. This method has 
been used by Griffiths and Warren (1952) for Co® and Na*™ gamma rays. 
West (1956), with much-improved apparatus and a detailed study of the over- 
all efficiency of the three-crystal spectrometer, was able to obtain absolute pair 
cross sections for the same gamma-ray cascades in both sodium iodide and 
anthracene. In the present work we have extended the energy range of the 
measurements of the relative pair production cross sections in sodium iodide 
by using Zn®, Na”, Co®, and RdTh sources of known strength in exactly the 
same geometry. Since some of these sources emit only single gamma rays above 
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1.02 Mev it was necessary to know the relative strengths of the sources in 
order to obtain relative pair production cross sections. The strength of the 
sources was independently measured to better than 5% with a sodium iodide 
scintillation counter employing a crystal 23 in. in diameter by 33 in. long whose 
efficiency as a function of gamma-ray energy had been studied in some detail 
in this energy range. (Details of this work will be published.) The ratio of pair 
production cross sections obtained for the gamma-ray energies given in 
brackets are: 

o (2.62) o(1.33)_ 9, o(t.l7) _ 


(1.33) sii wae o(1.11) 


If we take the absolute cross section for the 1.33-Mev gamma rays of Co*° 
from the work of West (1956) as (6.81+.17) X 10-°6 cm? then these ratios lead 
to the following absolute pair production cross sections in sodium iodide: 

(1.11) = (0.32+.03) X 10-°6 cm? 

(1.17) = (1.84+.12) X10-*6 cm? 

(2.62) = (116410) X10-*6 cm? 


4.22 


The errors of about 10% include that quoted by West for the 1.33-Mev 
gamma ray (about 2.5%) plus statistical and estimated systematic errors 
introduced by our ratios. The ratios obtained here for pair production by 
Co®, Na”, and RdTh gamma rays are in good agreement with the results 
collected by West (1956) and the Zn® result combined with the absolute values 


given by West for Co®’ gamma rays provides a measurement of the pair produc- 
tion cross section in sodium iodide at an energy about 90 kev above the pair 
production threshold. 
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POSSIBLE CHEMICAL REACTION OF ORTHOPOSITRONIUM WITH OXYGEN 


>: A... EW Pac. 


It has been shown previously (Paul 1958) that the effect of traces of oxygen 
in quenching orthopositronium in liquid argon is much greater than was 
observed by Deutsch (1953) when oxygen was introduced as the impurity in 
argon gas. 

In concluding the above-mentioned experiments (1958) two further measure- 
ments were made of the quenching of orthopositronium by oxygen at lower 
temperatures in liquid nitrogen. As before, the reciprocal of the r2 lifetime is 
plotted against oxygen concentration and a straight line can be fitted to the 
experimental points, which, however, are subject to rather large uncertainties. 
From the slope of the line the volume rate, A*/n, for the quenching process is 
obtained, A* being the increased annihilation rate due to the presence of n 
molecules/cm*® of Os. When a simple gas-kinetic picture of the quenching 
process may be assumed, one writes 


*. = oD (cm’ /sec) 


where o is the cross section for the quenching process and v the velocity of 
positronium. In liquids, however, it would seem plausible to make comparisons 
of volume rate rather than to use gas—kinetic ideas, where it is clear that they 
cannot be entirely appropriate. The new results together with the previous 
results are given in Table I. 


TABLE I 


Volume rate per mole- 
cule O2, 10- cm3/sec 
Medium Temperature, ° units Author 





Argon gas ~300 3.1 Calculated from 
Deutsch (1953) 


Nitrogen liquid 77.5 
Nitrogen liquid 64.7 


Argon liquid 87.5 g Paul (1958) 


Norte: The relative values of volume rates in liquid argon and liquid nitrogen are estimated to be accurate to 
within 20%. Absolute values are probably good to 20%. 


An electron exchange mechanism has been proposed by Ferrell (1958) to 
account for the high volume rate observed by Deutsch. It is probable, how- 
ever, that some additional mechanism must be invoked to explain the strong 
temperature dependence in the neighborhood of 80° K. 

A mechanism which would account for the more rapid quenching in liquid 
than in gaseous argon is the formation of the chemical compound (et e~)Os. 
That it is plausible that this positronium compound should be formed may be 
seen by considering its hydrogen analogue, hydroperoxo, the HO, radical. 
Hvdroperoxo, a stable but highly reactive compound, is readily formed from 
molecular oxygen and atomic hydrogen and has been detected with a mass 
spectrometer by Foner and Hudson (1953) and by Robertson (1954). Robertson 
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estimates furthermore that HO, formed through bimolecular association 
would dissociate in about 5X10-” second unless a collisional de-excitation 
should follow within that time. In three-body collisions, however, it would 
generally be possible for HO» to be formed in the ground state. If these argu- 
ments are applied to (e+ e~)Oz it is readily seen that a chemically stable com- 
pound would be formed much more abundantly, for example, in liquid argon 
with oxygen impurity than in the gas of identical composition, since the three- 











body reaction, 


(et e-)+0O.+X = (e+ e-)O.+X, 







would be quite rare in the gas, but probably every collision of (e+ e~) and Oy 
could be considered to be of this type in the liquid. In the gas phase posi- 
tronium would generally have undergone annihilation before a three-body 
collision would occur to form (e+ e~)O, in the ground state; and an excited 
state of (e+ e~)O2 would generally dissociate before annihilation took place, 
if one allows that the dissociation will not take longer, on the average, than 
a few times 10~" second. To account for the enhanced quenching in liquid 
argon, it must then be assumed that the ground state of (e+ e~)Osz is short- 
lived with respect to annihilation, having 7 ~ 10~-" second, so that the three- 
body chemical reaction is in effect competing directly with other modes of 
annihilation. This condition would surely be fulfilled, since the electron 
density at the positron must be at least as great in the compound as it is for 














free positrons in liquid oxygen. 

The total volume rate may thus be thought of as made up of a part due to 
annihilations which are consequent upon free collisions of oxygen and posi- 
tronium, and the rest a chemical reaction rate. The strong temperature depen- 
dence in the neighborhood of 80° K suggests fitting the data to an exponential 
function, and in fact the three experimental points will yield the constants in 










an expression of the type 






* 
z = A+Be**”, 


n 










It is tempting to identify the second term with the reaction rate, and obtain 
a value of ¢ ~ 0.09 ev, i.e., 2000 cal/mole. This quantity is a plausible activa- 
tion energy for the formation of (e+ e~)O2. The value of A obtained is then 
2.5+0.5X10-" cm/sec, less than but hardly different from Deutsch’s value 
for the quenching of orthopositronium by oxygen in argon gas. This last 
observation suggests that the term A arises from an electron exchange mecha- 
nism (Ferrell 1958) and that there is very little, if any, temperature depen- 
dence in the quenching arising from this mechanism. 
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NOISE TRAINS 
N. C. GERSON* AND W. H. Gossarpt 


A study of atmospheric noise in the frequency range 1-30 Mc/s was under- 
taken employing a dynamic spectrograph. The spectrograph included a set 
of R.F. comb filters whose outputs activated a moving strip recorder. The 
equipment thus provided an acceptance bandwidth of 500 kc/s with which 
noise bursts have been studied. During the investigation, it was found that 
some bursts appear to drift completely through the 500 kc/s frequency inter- 
val available. This note describes several of the preliminary results already 
noticed. Because of the frequency sweep of the atmospherics, the phenomena 
have been termed ‘‘sweepers.”’ 

When examined aurally with normal receiving equipment, the passage of a 
sweeper produces an audible swish, chirp, whistle, etc. The type of sound 
presented to the ear depended upon the bandwidth acceptance of the receiver 
and the rate of change of frequency with time (Af/A?’). 

The rate of sweeper occurrence in the range 24.0-24.5 Mc/s has been 
found to equal about 1250/hour between sunrise and sunset. The rate seems 
variable, and on other occasions the maximum rate exceeded 2300/hour and 
4200/hour, respectively. Generally it seems that sweepers at frequencies 
below 20 Mc/s are received mostly during darkness, while those between 20 
and 30 Mc/s, mainly during daylight. 

Sweepers already observed in the 500 Mc/s spectrograph cover frequency 
ranges greater than 75 kc/s. The total duration of an individual sweeper 
varies from 4 to 18 seconds and undoubtedly greater durations may yet be 
found. Sweepers may drift through 500 kc/s in less than 0.1 second. On an 
f-t graph, a typical sweeper begins as ¢ = constant, assumes a value of Af/A?? 
>0O, and terminates as f = constant. The time spent by the sweeper in each 
of these regions is variable. (Values of Af/At? = 5000 Mc/sec? have already 
been noted.) The great majority of sweepers drift from high to low frequencies 
although a few travel in the reverse direction. 

It has been found that sweepers appear in trains separated by a somewhat 
constant time interval. Trains composed of 50 to 100 individual sweepers in 
10 minutes have already been noted. Some trains appear to persist for hours. 
The resolution of the equipment has not permitted detection of significant 
changes in the dispersion of one sweeper relative to others in the same train. 

Sweepers probably originate as solar noise bursts although terrestrial elec- 
trical discharges or cosmic noise pulses must also be considered. The repetitive 
nature of sweepers in a train implies (if they are of solar origin) the existence 
of some relaxation process in the chromosphere. 


*Vice Chairman, Arctic Committee, United States National Committee for the International 


Geophysical Year. 
tDepartment of Defence. 
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If sweepers are of terrestrial origin, they may become trapped after radia- 
tion, by different magnetic force lines and be reflected successively between 
the northern and southern hemispheres. In this case, they become an extension 
of whistlers, except that they occur at much higher frequencies. The periodic 
radiation of energy may occur from bunched electrons (either in the earth's 
atmosphere or in interplanetary space) moving in crossed electric and magnetic 
fields. 


RECEIVED May 1, 1959. 
PRIVATE CONTRIBUTION. 


AUTOCORRELATION OF QUANTUM-MECHANICAL WAVE FUNCTIONS 






RICHARD BOURRET 








I. CORRELATION FIELD THEORIES* 


In recent years autocorrelation functions have received a great deal of 
attention in such disparate studies as statistics, the theory of turbulence, 
electrical filter analysis, optics, and acoustics. One of the early proposed 
applications was by Richardson (1926), who suggested the use of the auto- 
correlation of diffusate concentration instead of concentration itself in the 
differential equation governing turbulent diffusion. In an equation for turbulent 
diffusion recently advanced by the author autocorrelation functions play a 
central role, but in this case they are autocorrelation functions of turbulent 
velocity fluctuations rather than of the concentration of the diffusate. These 
velocity autocorrelations are the basic variables in the modern theory of 
turbulence (e.g. Batchelor 1956). In acoustics the autocorrelation of the 
pressure field has been usefully employed where noise contamination has 
been present (e.g. Lyon 1950). In the optical theory of partially coherent 
radiation autocorrelation functions of the electrical and magnetic fields have 
been used by Wolf (1954, 1956) and later by the author (Bourret 1959). 
In these studies the autocorrelation functions are shown to satisfy “‘field’’ 
equations formally identical with the Maxwell equations. 

In view of this variety of applications of correlation functions to physical 
theory, it seems worth while to consider the possibility of a quantum-me- 
chanical autocorrelation function constructed from the state function of a 
system. 


























II. DEFINITION OF AUTOCORRELATION FIELD IN QUANTUM 
MECHANICS 
We begin with a definition of the autocorrelation field of a system in an 
eigenstate of the energy. Designating the associated wave function by Wz(x) 
we define the complex autocorrelation field by 







*This material was done under contract Nonr 840(oi) and is contribution No. 235 from the 
Marine Labaratory, the University of Miami. 
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NOTES 









(1) x(n) = f Vale) ¥e(e+n) dx. 


Some obvious properties of this function are: 






(2) oe(—7) = on(n) 
oz(n) - 0 nto 










and 
(4) on2(0) = 1 


on account of the assumed normalization of the wave function. An interesting 
consequence of our definition is that shifts in the co-ordinate origin of the 
wave function will leave the autocorrelation function unchanged. Since both 
quantum mechanics and relativity forfeit classical information in exchange 
for applicability to non-classical problems, and since general relativity, in 
particular, forfeits ‘‘co-ordinate knowledge’”’ altogether, it may be hoped that 
the information lost in the transition to the autocorrelation function may be 
physically inessential. 

More knowledge is lost, however, in this transition than mere knowledge 
of the co-ordinate origin as can be seen from the following: from the usual 
relation between the wave function and its momentum transform 














(5) xe(p) = h J vec FI dee 





we obtain by inversion and substitution into equation (1) 





(6) oz(n) = f |xe(p)|* e?""dp. 





In the absence of a magnetic field the wave function may be chosen real, and 
so equation (5) implies that 

si & 

(7) xe(—p) = xe(P). 


We can, therefore, invert (6) and express it as 







(8) xe(p)x2(—p) = 0 f dx (nye Man. 





Because of our specialization to energy eigenstates ~z(x) can be recovered 
from ¢z(n) providing only that the Hamiltonian be reflection invariant, for 
in that case the wave function must have even or odd parity (Schiff 1949) and 
hence, by equation (5), so must xz(p). Equation (8) then shows that z(x) 
is determined to within a constant phase factor. 









III. OTHER PROPERTIES OF THE AUTOCORRELATION FIELD 


From equation (6) it is clear that ¢g(n) is simply the expectation value of 
the dynamical variable e*”/". This can also be seen by expressing the defining 


relation, equation (1), in the form 
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(9) $e(n) = J Va xe?" Ya(x)dx = (e”"")p 


and regarding p as an operator. From this last relation we see that ¢z(n) can 
be employed to evaluate the expectation value of any function of the momen- 
tum by the operation 


(10) Gierre = 12 2) dala) lne 


Because the preceding formula can readily be used to obtain the moments of 
the momentum distribution we see that ¢z(n) corresponds to what is called 
the moment-generating function in statistics. 

An interesting feature of the interpretation of ¢g(n) as an expectation value 
of a dynamical operator is that it enables us to construct its exact classical 
analogue. If P(p) is the classical momentum distribution, we can for the classi- 
cal ‘‘autocorrelation function equivalent” by forming the function 


(11) x(n) = fe™"P(p)dp. 
We may, alternatively, construct it as a time average by the formula 
i; 
(12) oz(n) = T eth 
0 


where T has been suitably defined. 

The preceding definitions of quantum-mechanical autocorrelation function 
for systems in an energy eigenstate can readily be generalized to general state 
functions. Consider the general time-dependent state function 


(13) v(x, t) = A} fel valeye dE 


where |c(£)|? is the probability density of the energy states. We generalize 
our definition simply by including an autocorrelation over time, i.e. 


Si F + 
(14) o(n, 7) = f at fax wv (x, t)W(x+n, t+7). 


Substitution of (13) into (14) yields 


(15) (nt) = J c(E) "a(n *"GE. 


This can be replaced by a sum when discrete energy levels are to be used. 
The case in which r = 0 takes the specially simple form of a weighted average 
of the ¢g(n). The classical analogue of (15) is obtained by replacing E with 
the Hamiltonian for the system and integrating over phase space. The weight- 
ing function |c(£)|? is then replaced by the normalized Gibbsian expression 


(16) exp[— H(p, q)/6] 
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to give the canonical autocorrelation function analogue. 6, as usual, is the 
modulus of the distribution and is defined by 


(17) 6= kT. 
IV. EXAMPLE OF THE AUTOCORRELATION FIELD, QUANTUM AND 
CLASSICAL 
The harmonic oscillator provides a convenient illustrative example of the 
autocorrelation field. Using the abbreviations 
(18) B? = mw/h 
and 
w = classical angular frequency of oscillator, 
we write the wave functions corresponding to the energy 
(19) Ey = hw (n+3) 
as 
(20) Yn(x) = (2"n!)3(8?/m)* eH, (Bx), 
where the H, are Hermite polynomials. Using these functions in the defining 


equation (1), and employing an expansion theorem for Hermite functions 
(Morse and Feshbach 1953) we obtain for the autocorrelation field of the 


nth eigenstate 


\ 1 —B2n2 * 4 
(21) oa(n) = 3 €* “Ly(6'"/2) 


where L, is the Laguerre polynomial of order n. It is interesting to note that 
these functions are the axially symmetric solutions of the two-dimensional 


Schrédinger equation 

‘ 2 
(22) voto a 7) ¢=0 

h 2 
with the radial distance 7 replaced by 9/2, and with the normalization appro- 
priate to autocorrelation functions. This relation to the two-dimensional 
harmonic oscillator gives some additional insight into the quality of the in- 
formation lost in the transition from the wave to the autocorrelation function. 

The classical autocorrelation function equivalent may be derived from 

expression (12), taking 7 to the oscillatory period. For the harmonic oscillation 
described by 
(23) x = xo sin w (t+6) 


and 
p = mxww cos w (t+64) 


where x» and 6 are arbitrary, equation (12) readily yields 


6 3 
(24) oz(n) = ee) *| 
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where we have used 
(25) E = 3mwx, 
and Jo is the Bessell function of order zero. 
The canonical autocorrelation function for a one-dimensional assembly of 
harmonic oscillators with a Boltzmann distribution in energy (the Einstein 


solid) is constructed from equation (15), rewritten as a discrete sum, with the 
weighting factors 


lI 

7 
& 
= 

= 


(26) ree where Zz 


Using the notations 


] 1 ir hw 
27 ee —_— . = < ciate 
(27) O(n) Ot and A(r) = tanh 26(7) 
the calculation yields 
ia _ _sinh(hw/8) bi ma | 
me) O(n 7) = sinh {hw/0(r)} PL ahACr) I 


The special case with r = 0 and with high temperature will be of interest. 
It is simply 

—mOn2 /2h2 
(29) #(n:0) = ener 
where 


hw 
20 at. 


The corresponding classical canonical function is obtained from equation 
(15) by using Gibbsian weights and replacing the energy integral with an 
integral over phase space. We have: 


0) nr) = 2" f faxap enone 5] 2 /imliD, 5) | 


where Z now has the classical significance 


(31) Z= J facap foe 


Substituting the harmonic oscillator Hamiltonian 





‘ p mw 2 
9 4 \= adem Aim 
(32) H(p, *<) 2 + 9 x 





into equation (30) we obtain 


_ Or) E43 
(33) o(n, 7) = 9 xP OR : 


For r = 0 this becomes identical with the high-temperature result for quan- 
tum-mechanical assemblies. 









NOTES 


V. CONCLUSION 


From the foregoing it is clear that the autocorrelation function for a 
quantum-mechanical system is intermediate between the state function and 
the partition function in the amount of physical information which it carries. 
If, therefore, further investigations should provide us with a means of calculat- 
ing it directly (i.e., not via the wave function nor, classically, via the solutions 
of the equations of motion) it might prove a useful tool when the complete 
wave-mechanical description is inaccessible and the partition function in- 
adequate. 
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ANHARMONIC FORCES AND THE EINSTEIN MODEL OF A CRYSTAL* 


J. W. LEEcHT 


At temperatures near the melting point the specific heats of crystals are 
sometimes found to exceed the classical values predicted by the harmonic 
model (see for example Blackman 1955; Berg and Morrison 1957). Such anoma- 
lous values have been attributed to the action of anharmonic forces. 

The mathematical problem of allowing for anharmonic forces is a formidable 
one and it is not surprising to find their effect calculated in an approximate 
manner using the simple Einstein model of a solid (see Henkel 1955; Zucker 
1958). 

The validity of such an extended Einstein model could be estimated by 
comparing its predictions with experimental observations. It is doubtful, 
however, whether existing experimental data are sufficiently accurate in view 
of alternative explanations in terms of lattice defect formations (see for 
example Carpenter 1953; MacDonald 1953). The predictions could also be 
compared with those of a fuller mathematical analysis based on a more 


*Issued as N.R.C. No. 5299. 
tNow at the Department of Physics, Queen Mary College, University of London, London, 


England. 
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obviously justifiable model. There have been such theoretical approaches 
(see for example Born and Brody 1921; Schrédinger 1922; Peierls 1929; Born 
1933; Hooton 1955; Ludwig 1958; Stern 1958) but the conclusions, aside 
from general predictions concerning temperature dependence, are not yet in a 
conveniently usable form. In the absence of direct checks on the validity of 
the anharmonic Einstein model it is useful to make an indirect one by con- 
sidering the one-dimensional case, which, though unrealistic, has the merit 
of being mathematically manageable. 

An exact classical statistical mechanical analysis of the anharmonic linear 
chain model, considering only nearest-neighbor interaction, has been made by 
Damkohler (1935), by Dugdale and MacDonald (1954), by MacDonald and 
Roy (1955), and by Leibfried (1955). Their values for the excess specific heat 
at constant length may be expressed as: 


(1) AC, = ———— (As"— 204) +0(T*) = CT(As*—24) 


where D is the mutual P.E. of two atoms at their equilibrium separation ro, 
and the potential energy at separation r = (ro+x) relative to that at ro is 
given by 


V = D{ (ax)*+d3(ax)"+Aa(ax)*+ .. «4. 


Damk3ohler’s result differs from the others in the coefficient of the second- 
order term, but this is immaterial in the present analysis. 

In the Einstein model approximation each atom of the chain moves in the 
field of its two nearest neighbors considered frozen at their equilibrium 
positions. Effects due to odd-order terms of the potential series cancel out 
and a short analysis shows that to first order in 7: 


AC, = CT(—2)4). 


In the specific forms normally assumed for the interatomic potential 
(Morse and Lennard Jones—Mie types) \4 is positive and comparable in magni- 
tude to \;*. If we regard (1) as the standard result for the one-dimensional 
problem, it follows that the cancelling out of the odd-order terms in the 
analysis of the anharmonic Einstein model leads to serious error in AC,. 

There is no guarantee that one-dimensional considerations have any bearing 
on the three-dimensional problem, but some sort of parallel would be ex- 
pected. On this basis it appears unlikely that the anharmonic Einstein model 
will provide reasonable values for the specific heat. This conclusion is at 
variance with the apparent agreement with experiment of the conclusions of 
Henkel and Zucker. Domb (private communication) has suggested that one 
would expect a ‘mean’ field approximation to be reasonable in three dimensions 
since the number of nearest neighbors is relatively large. In support he points 
out that in the corresponding approximation for the Ising model the results 
are quite wrong in one dimension but reasonable in three dimensions. If this 
suggestion can be justified in detail it points to the need for caution in extra- 
polating from conclusions based on one-dimensional analysis. Meanwhile 
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there still seems to be doubt whether the observed excess specific heats can 


be attributed to anharmonic forces. 
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THERMAL CONDUCTIVITY OF SOLID ARGON AT 80° K 


D. J. LAwrence, A. T. STEWART, AND E. W. GUPTILL 


We have measured the thermal conductivity of solid argon by a simple 
technique which is suitable for about 10% accuracy. Figure 1 shows a schematic 
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Fic. 1. Schematic diagram of the experimental arrangement. 
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diagram of the apparatus. Commercial argon flows at a constant rate through 
a wet-type gas meter and a drying tube into a 2.5-cm diameter precision bore 
glass tube immersed about 26 cm in liquid air. The argon freezes in a clear 
polycrystalline layer of uniform thickness on the inside of the glass tube. The 
heat of sublimation is conducted away through the solid argon and glass layers. 

The heat flow, Q, in a radial direction through a uniform cylindrical shell of 
thermal conductivity, k, and length, h, is given by 


Q = —k2arh dT/dr 
where d7°/dr is the temperature gradient at a distance r from the axis. Integra- 
tion of this equation through the glass tube and the solid argon shell and 
elimination of the interface temperature gives 
2ahAT_ 1 
Q ka 
where AT is the temperature difference between the solid argon inner surface 
and the bath. The second term on the right is constant involving only the 
dimensions and conductivity of the glass, so the solid argon conductivity, ka, 
can be obtained from the slope of a graph of 2mhAT/Q versus In (r,/r). The 
heat, Q, which flows through the solid argon is obtained mostly from the heat 
of sublimation of solid argon. We used a value of 189 joules/g, which is 
obtained from the Clausius—Clapeyron equation using the vapor pressure 
measurements of Clark et al. (1951). The remainder of Q (up to 10%) is supplied 
by gas cooling. The amount of heat is calculated from the thermal data of 
Hilsenrath et al. (1955) and the measured temperature change of the gas. The 
initial temperature was measured just above the solid argon top edge by a 
thermocouple suspended in the flowing gas. The final temperature of the gas 
was that of the inner surface of the solid argon, which was obtained by vapor 
pressure measurement. A plot of solid argon vapor pressure for a typical run is 
shown in Fig. 2. The flow was stopped three times in this run to enable the 
slowly changing bath temperature to be measured. The pressure fluctuations 
in this run are due to occasional sticking of a bath stirrer. These readings 
were used to calculate the values of AT shown in Fig. 3. The corresponding 
values of In (7;/r) were calculated from the gas meter volume measurements and 
the solid argon density (Dobbs and Jones 1957). From each of the three 
sections of Fig. 3 we obtained a value for the slope, and from the weighted 
average the value of thermal conductivity which is shown. 

The error shown is due only to internal consistency in vapor pressure 
measurements. The systematic uncertainties in the experiment are of the 
order of 15%. The final average value from seven runs for the thermal con- 
ductivity of argon at 80° K is 


In (n/t In (72/11) 
G 


k = 2.6+0.4 mw/cm deg. 


This is in reasonable agreement with White and Woods (1958) value, 
k = 3.140.3 mw/cm deg (taken from their graph), who used quite different 
techniques of formation of the solid and of measurement. 
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ARGON VAPOR PRESSURE (CM HG) 
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Fic. 2. Typical argon vapor pressure data for one run. 


Fic. 3. Plot of log (r:/r) is a function of AT. The argon conductivity ka is obtained from 
the slope. 
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THE DECAY OF Mg” 
R. S. StorEyY* AND K. G. MCNEILL 


Magnesium-23 decays by positron emission with a half-life of 12 seconds. 
This decay has been studied by a number of authors (Hunt 1954; Boley 1951) 
and the energy of the beta transition to the ground state of Na* has been found 
to be 2.95 Mev. The existence of an energy level in Na® at 0.439 Mev above 
the ground state has been shown by work on the decay of Ne® (Penning 1957) 
and on coulomb excitation of Na** (Temmer 1956). A positron decay of Mg” 
to this level is energetically possible, and moreover is allowed on Gamow-Teller 
selection rules ((3/2)+—>(5/2)+). Such a transition should be followed by a 
0.44-Mev gamma ray. In 1939 no gamma ray was observed during experiments 
on this nuclide (White 1939), but modern techniques are much more sensitive 
and discriminatory than the absorption methods used by White et a/. In view 
of the recent interest in this isobar (Paul and Montague 1958), a search for the 
0.44-Mev gamma ray has been made using scintillation spectrometry. 

Magnesium, 99.8% pure, in the form of a cylindrical lump 1}-in. diameter 
by 2 in. long, was irradiated in the bremsstrahlung beam of the 23-Mev betatron 
of the Ontario Cancer Institute, Toronto. After irradiation for 15 seconds, the 
magnesium was quickly transferred to a counting position. Here gamma rays 
from the magnesium could be detected by two 2 in. X 1} in. Nal crystals 
mounted on DuMont 6364 photomultipliers. The two detectors made an angle 
of 80° to each other, the Mg source being at the point of intersection of the 
axes of the photomultipliers. Scattering from one crystal to the other was 
prevented by lead shielding. One of the detectors fed a 100-channel pulse height 
analyzer (100 PHA). The other fed a single-channel PHA, set to pass pulses of 
0.48 to 0.54 Mev. Output pulses from this gated the 100 PHA. Thus into the 
100 PHA would pass only those gamma rays in coincidence with positron 
annihilation quanta. Moreover, the angular separation of the detectors pre- 
vented any true coincidences between the annihilation quanta. The resolving 
time of the coincidence circuit was 1.6 microseconds. 
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Counting lasted 15 seconds. After a lapse of a further 15 seconds, a back- 
ground of 15 seconds was subtracted. The sample was then reirradiated. This 
cycle was repeated until a sufficient number of counts was accumulated. The 
removal of background after each count eliminated the effect of build-up of 
background due to 15-hour Na™ made by the reaction Mg**(yp)Na*. The 
results of one such run, as displayed on the 100 PHA, are shown in the figure. 


Counting Rate 


0 
o 


aaa __- 





O-3 o4 os 6 
Energy in Mev 


Fic. 1. The spectrum of gamma rays in coincidence with annihilation quanta. The 0.51-Mev 
peak is due to random coincidences. The dotted line indicates the spectrum obtained from C", 
a pure positron emitter, the areas of the 0.51-Mev photopeaks being the same in both cases. 
The energy scale is calibrated using Na®? and Cs'3?7 sources. 


In the coincidence spectra displayed on the 100 PHA, there were two photo- 
peaks apparent. One, at 0.51 Mev, was due to random coincidences with 
annihilation quanta, whereas the other, at 0.44 Mev, could not be due to this 
cause—this fact was checked using C", a pure positron emitter. We attribute 
it to the decay of Mg** by the emission of a 2.51-Mev positron followed by a 
0.44-Mev gamma ray from the first-excited state of Na”. 

An estimate of the relative intensity of the 2.51-Mev positron decay compared 
with that of the 2.95-Mev decay was made, using the results of the coincidence 
measurement and the fact that without the coincidence circuit operating, only 
a slight distortion of the low energy tail of the 0.51-Mev peak was seen at 
0.44 Mev. The result is that 0.065+.025 of the disintegrations are via the 
excited state of Na”. The inaccuracy quoted is largely due to difficulties of 
beam control caused by the cycling procedure. 

Magnesium-23 has a half-life of 12 seconds. Using the experimental result 
that 6.5% of the disintegrations are to the excited state of Na*, the partial 
half-lives for the ground states — excited state transitions are 12.7 seconds and 
185 seconds, with logy ft values of 3.7 and 4.5 respectively. The first value is 
characteristic of a super allowed transition, caused by the fact that Mg** and 
Na*™ are mirror nuclei of each other in their ground states. The value of 4.5 is 
consistent with the expected allowed character of a 2.51-Mev positron 


transition. 
To summarize, the existence of an allowed transition from the ground state 


of Mg” to the 0.44-Mev excited state of Na”* has been demonstrated. 
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EXCHANGE REACTION OF 90-MEV PROTONS ON DEUTERONS 
D. R. FirtTH 











INTRODUCTION 


Exchange forces are demonstrated by the appearance of high-energy protons 
when hydrogen is bombarded by neutrons of energies much greater than the 
potential well (Fermi 1950). 

In like manner one expects fast neutrons from the reaction of fast protons 
on deuterons. Indeed, predictions for energy and angular distributions of such 
neutrons are implied in the calculations of Gluckstern and Bethe (1951) and 
Chew (1951) on the mirror reaction at 90 Mev. Of especial interest is their 
prediction of a narrow peak only 1 Mev wide at the maximum neutron 
energy. This maximum energy decreases as the scattering angle increases. 

However, the spread of at least a few million electron volts commonly 
found in the proton beams, as well as the relatively poor resolution of neutron 
detectors at 90 Mev, means that one cannot accurately check this process by 
looking at the neutrons. This was demonstrated in the work of Hofmann and 
Strauch (1953), who attributed the shape of their neutron spectrum to a 
spread of 25 Mev in proton beam energy. 

On the other hand, one may look at the total momentum and total 
energy of each pair of remaining protons, and from conservation laws deduce 
the energy and angle of the emitted neutron. The shape of the neutron spec- 
trum is not expected to vary much from, say, 85 Mev to 90 Mev, so that 
the spectrum of the total pair energy (which is the reverse of the neutron 
spectrum taken to the origin) will not be blurred by a spread in beam energy. 
If the incident energy is somewhat uncertain there will be only a small error 
in the calculated angle of the neutron. 

The probability of this reaction, however, will be strongly influenced by 
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PLATE I 





Fic. 1. As an example of the analysis, event la shows two protons of energy 5.6 Mev and 
5.8 Mev. Taking the average incident proton energy as 87 Mev, one finds the neutron emitted 
at 29°. From Gluckstern and Bethe, for an 87-Mev proton and the neutron off at 30°, the 
most likely energy for the proton pair is 11.4 Mev, which agrees well with (5.6 plus 5.8) Mev. 
In its center of mass this proton pair has an energy of 1.4 Mev. 
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the final state interaction between the slow particles (Watson 1952; Migdal 
1955). One can imagine, for example, a primary mechanism of charge exchange 
between the incident proton and the neutron, leaving two protons with low 
relative energy, which will therefore have / = 0 and S = 0. But in the 
original deuteron these two nucleons have / = 0 and S = 1, so the charge 
exchange interaction must also change the spin state. The final state inter- 
action is also largely responsible for the sharpness of the neutron spectrum; 
final states are favored which have the two protons within p-p virtual 
resonance. 


EXPERIMENTAL 


The aim of this note is to show the feasibility of investigating this reaction 
by means of a diffusion cloud chamber, which was filled to 19 atmospheres 
with deuterium. 

Pulses of protons at 90 Mev from the external beam of the McGill synchro- 
cyclotron were sent through the cloud chamber and photographs were taken 
120 milliseconds later. It was possible to use relatively intense proton beams, 
since the slow protons from the reaction give tracks more than 20 times as 
dense as 90-Mev protons, and therefore stand out among the mass of tracks. 

The photographs (Fig. la and 1b) show two examples of the charge 
exchange. Preliminary results from several hundred photographs indicate that 
about one out of every six events is due to charge exchange. 

In both figures certain bare patches appear in the beam tracks. These are 
caused by the passage, some seconds before the proton beam, of an ionizing 
particle which removes vapor from a considerable region. Another feature 
typical of a diffusion chamber is the uneven supersaturation throughout 
the sensitive region, which can cause a track to be less dense at its end than 
its beginning (see event on Fig. 1b). 

It is concluded that the neutron—proton exchange force is spin dependent 
and can be investigated in some detail in this way. 
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MOMENTUM DISTRIBUTION OF ANNIHILATING 
ELECTRON-POSITRON PAIRS IN LiF 


R. H. Marcu AnD A. T. STEWART 


We have determined the momentum distribution of electron—positron pairs 
annihilating in LiF in order to resolve the large discrepancy (see Figs. 1(B) 
and 1(C)) which exists between the results of Lang and DeBenedetti (1957) 
and Millett and Castillo-Bahena (1957). Our results agree fairly well with 


Lang’s. 





Fic. 1. Momentum distribution of electron—positron pairs annihilating in LiF: (A) this 
experiment, (B) the distribution obtained by us from Lang’s smoothed data (1957), (C) as 


given by Millett (1957). 


Can. J. Phys. Vol. 37 (1959) 
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Before presenting the various data it is necessary to discuss briefly' the 
experimental arrangements. The technique is the usual one of measuring the 
angle between annihilation photons to determine their momentum and hence 
the momentum of the annihilating pair. Lang’s apparatus and ours detected 
annihilation radiation in two scintillation counters mounted behind slits in 
Pb shields. The source, midway between the counters, was arranged to move 
at right angles to the plane of the slits. Coincidences were observed as a 
function of the transverse displacement of the source. The rate of coincidences 
between the two counters is proportional to the number of events which 
have a given z component of momentum, p,, determined by the displacement 
of the source (taken to be along the z axis). 





° 2 4 6 8 10 12 14 px 10 


Fic. 2. Comparison of the momentum distribution of electron—positron pairs annihilating 
LiF and NaF. @ LiF, O NaF (Stewart and Pope). 


Millett’s apparatus defined the angle between the annihilation photons by 
accepting in one counter a narrow pencil of radiation through a small circular 
aperture, and in the other a cone of radiation defined by an annular aperture 
in Pb shielding. The angle between accepted photons was varied by inserting 
annular apertures of different sizes, the relative efficiencies of which must be 
known. The rate of coincidences between the two counters in this experimental 
arrangement is proportional to the number of events that have a given radial 
component of momentum, determined by the angle of the cone. (The polar 
axis is taken to be the line joining the counters.) From both experimental 
arrangements the momentum distribution, N(p), may be obtained by an 
inversion process. 

Momentum distributions obtained from an inversion of our data, from an 
inversion of Lang’s data, and by Millett are shown in Fig. 1. It is seen that 


1A full discussion of these experimental arrangements can be found in the references and 
A. T. Stewart, Can. J. Physics, 35, 168 (1957). 
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Lang’s data gives a slightly higher momentum for the peak of the distribution 
than ours. There is, however, a large discrepancy between both these distri- 
butions and that obtained by Millett and Castillo-Bahena. It is possible that 
such a discrepancy could occur if the relative efficiencies of the various annular 
apertures used by Millett were not accurately known.? This problem is not 
encountered in the parallel slit arrangement because one slit is used for all 


angles. 

The result obtained for LiF is compared with a similar result for NaF? in 
Fig. 2. The two are similar thus confirming the commonly held hypothesis 
that most annihilations occur with the electrons of the negative ion. 


NOTE ADDED IN PROOF.—Recent data obtained by Dr. Millett is in reason- 
able agreement with results presented here. We are indebted to him for his 
courtesy in sending us the new results in advance of publication. 


We are grateful to J. M. Robson for the assistance given us in connection 
with the shipments of Cu®™ and thanks are due to Atomic Energy of Canada 
Limited for the loan of equipment. 
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MILLETT, W. E. and CastTILLo-BAHENA, R. 1957. Phys. Rev. 108, 257. 


RECEIVED JUNE 2, 1959. 

DEPARTMENT OF PuHysICcs, 
DALHOUSIE UNIVERSITY, 

HALiFax, N.S. 


2It may be noted that humps occur at the same angle, about 6 milliradians in Millett’s data 


for all materials. 
3We are indebted to Stewart and Pope for these results in advance of publication. 
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